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Abstract. We propose a definition of center of mass for asymptot- 
ically flat manifolds satisfying Regge-Teitelboim condition at inflnity. 
This definition has a coordinate- free expression and natural properties. 
Furthermore, we prove that our definition is consistent both with the 
one proposed by Corvino and Schoen and another by Huisken and Yau. 
The main tool is a new density theorem for data satisfying the Regge- 
Teitelboim condition. 



1. Introduction 



A three-manifold M with a Riemannian metric g and a two-tensor K is 
cahed a vacuum initial data set (M, g, K) if g and K satisfy the constraint 
equations 

Rg-\K\l + H^ = 0, 
dwg{K) -dH = 0, (1.1) 

where Rg is the scalar curvature of M and H = Trg{K) = g^^Kij. We say 
{M,g,K) is asymptotically flat (AF) if it is a vacuum initial data set and 
there exists a coordinate {x} outside a compact set, say B^^, in M such 
that 

gij{x) = 5ij + hij{x), hij = 0(|x|~^) Kij{x) = 0{\x\~^) 

rnjM^) = 0(1x1-2) ^ 0(1x1-3) 

9ij,kii^) = 0(|x|-3) K,j,ki{x) = 0(|x|-^), (1.2) 

and similarly for higher derivatives. For AF manifolds, the ADM mass m is 
defined by 

{9ij,i - 9ii,j)^i dag, 

where {|x| = r} is the Euclidean sphere, Vg is the unit outward normal vector 
field with respect to the metric g, and dag is the volume form induced from 
iM,g). 
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For the case in which we are interested, we require asymptotic symmetry 
on {M,g,K). We say {M,g,K) is asymptoticahy flat satisfying Regge- 
Teitelboim condition (AF-RT) if (M, g, K) is AF and g, K satisfy these 
asymptoticahy even/odd conditions 

5°f'^(x) = 0(|x|-2) i^--(x) = 0(1x1-3) 

{g^'')^,{x) = 0(|x|~3) {Kfr),,ix) = Oi\x\-'), (1.3) 

and on higher derivatives, where f°'^'^{x) = f{x) — f{—x) and f'^'"^^{x) = 
fix) + /(-x), [RTTi] . Notice that f"'^'^ and r"^"" are only defined outside 
B/Jq in which the coordinate is defined. It is proved by Corvino and Schoen 
in [CS06] that AF-RT manifolds form a dense subset of AF manifolds in 
some suitable weighted Sobolev space. 

For (M, 0, K) satisfying AF-RT, we propose an intrinsic definition of cen- 
ter of masqj suggested by Richard Schoen, 

= le^.^-l^X,., ^^'^ - \Ra9^Mc^)^ld^3, a = 1,2,3 (1.4) 

where Rij is the Ricci curvature of M and Y(^) = (IxpJ'^* — 2x'^x*)^ is a 
Euclidean conformal Killing vector field. If {M,g,K) is AF-RT, the limit 
converges. Notice that if we replace Ug by the normal vector field with 
respect to the Euclidean metric and dug by the volume form with respect 
the induced metric from the Euclidean space, the limit is the same. Also 
notice that the above expression is not defined when the ADM mass m is 
zero. However, the center of mass can not be well-defined when m is zero 
because a basic version of positive mass theorem shows that M is actually 
the Euchdean space [SY79] [SYST] . 

This intrinsic definition is motivated by a similar expression of the ADM 
mass when Y(q,) in the above definition (jl.4p is replaced by {— 2x*^}, the 
radial direction Euclidean conformal Killing vector field, 

If 1 

m = lim / (Rii Raqij]{—2x^)vldaa- (1-5) 

167r r->oo _/|^|^/ 2 ^^^'^ 3 ^ ^ 

The Euclidean conformal Killing vector fields {— 2x*^} and Y(„) generate 
dilation and translation at infinity. A detailed discussion about the cor- 
respondence between these two vector fields and the mass and center of 
mass can be found in |C W08j . Another motivation comes from the spatial 
Schwarzschild metric 

= 1 + 17, 1 = 1 + ^ + ^13^ + 0(|x|-3). 

\ 2\x—p\J \x\ \x\'^ 2\x\-^ 

^The intrinsic definition can be generalized to n-dimensional manifolds if we replace 
the factor by a suitable constant depending on n and LOn, where ujn is the volume of 

the unit ball in R". Hence the following arguments work more generally for n-dimensional 
manifolds, but for simplicity, we assume n = 3. 
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If we replace the metric in ()1.4p by g'^ , Cj is precisely the vector p which 
indicates the center of the manifold. It is worth mentioning that although 
in the Schwarzschild case p is not a point in the manifold, one can develop 
polar coordinates using concentric spheres centered at p. 

Other definitions of center of mass have been proposed. Huisken and 
Yau |HY96j define the center of mass for (M, g, K) which are spherically 
asymptotically flat (SAF), i.e. {M,g,K) is AF and 

gijix) = (i + j^) +Pij^ M^) = o{\x\-^),dyij{x) = odxr^-hi)^ 

(1.6) 

They first prove the existence and uniqueness of the constant mean curvature 
foliation {Mr} for SAF manifolds, where the mean curvature of Mr is ^ — 
^ + 0(r~^). Then they show that the approximate centers converge to the 
center of mass Chy defined as follows: 

r. , dao 

C^y=lim^f— , a = 1,2,3, (1.7) 

where z is the position vector of Mr in M. This definition is also motivated 
by the spatial Schwarzschild manifold in which the constant mean curvature 
foliation is {\x — p\ = r}. Another definition by Corvino and Schoen in 
[CS06| is defined for AF-RT manifolds by 



\x\-- 



{hiaK - hiiv")dag , (1. 



where we recall that hij = gij — 5ij. One application of this definition in 
|CS06| is the gluing theorem which allows Corvino and Schoen to approxi- 
mate AF manifolds by solutions which agree with the original data inside a 
given region and are identical to a suitable Kerr solution outside a compact 
set. Their definition is a local-coordinate expression which is convenient for 
calculation purposes, but obscures interesting physics and geometry. 

The main purpose of this article is to prove that the intrinsic definition 
p.4p is equivalent to the Corvino-Schoen definition (|1.8|) . Moreover, for SAF 
manifolds in which the unique foliation of constant mean curvature surfaces 
exists, our intrinsic definition is equal to the Huisken- Yau definition. In 
other words, the intrinsic definition is a coordinate- free expression of the 
Corvino-Schoen definition and it generalizes the Huisken- Yau definition. 

Theorem 1. Assume {M,g,K) is AF-RT / fO) . Then 

Ci = Ccs- 
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Theorem 2. Assume (M, g, K) is SAF / fiT^) . Then 

Ci = Chy- 



We would like to make a note that Corvino and Wu jCWOSj recently have 
a result about the equivalence of these definitions under the assumption that 
the metric g is conformally flat at infinity with vanishing scalar curvature. In 
that special case, they are able to derive some explicit estimates. However, 
it seems that their approach cannot be generalized to AF-RT metrics. 

This article is organized as follows. In section 3, we prove a density 
theorem (Theorem 13. 2p for AF-RT manifolds. The theorem is crucial for 
most of the arguments in the article and may be of independent interest. In 
section 4, we discuss properties of the intrinsic definition. Theorem [1] and 
Theorem [2] are proved in sections 5 and 6 respectively. 
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3. The Density Theorem 

Let (M, g, K) be a vacuum initial data set. We introduce the momentum 
tensor 

^ij = K'^ -Trg{K)g'^ . 
The constraint equations (II. ip then take the form 

divg(^) = 0, (3.1) 

and we define 

In the case that {M,g,TT) is AF-RT, Corvino and Schoen jCSOGj prove that 
AF manifolds can be approximated by AF-RT manifolds in some weighted 
Sobolev space. Moreover, instead of requiring smooth solutions to the con- 
straint equations (j3.ip . their theorem works for solutions with weak regu- 
larity. Before we state the theorem, we need the following definitions. 
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Definition 3.1 (Linear and Angular Momentum). The linear momentum 
(Pi, P2, P3) and the angular momentum {Ji, J2, J3) are defined as follows. 

1 f 

Pi = — lim / TTiii^i da 



11111 I "11'^ a "''-'gi 



' |a;|=r 

where give the rotation fields in , for example = x^-^ ~ ^^gfi^- 

Remark. The linear momentum is well-defined for AF manifolds and the 
angular momentum is well-defined for AF-RT manifolds. 



Definition 3.2 (Weighted Sobolev Spaces). For a non-negative integer k, 

-5 



a non-negative real number p, and a real number S, we say f G w'^f{M) if 



W'lfiM) 



|a|<fe 

where a is a multi-index and p is a continuous function with p = \x\ on 

M\Br,. 

When p = CO, 

ll/llTy^-(M) = E esssup|Z)"/|pH+<5. 

^ ' rz, M 

\a\<k 

Definition 3.3 (Harmonic Asymptotics). (M, (7, vr) is said to have harmonic 
asymptotics if (M, g, n) is AF and 

g = u^6, TT = u^{CsX) (3.2) 

outside a compact set for some u,X tending to 1,0 respectively, where for 
any metric g, CgX is the operator associated to the Lie derivative Lxg 
defined by 

CgX = Lxg - divg{X)g. 

By the constraint equations ()3.ip , u and X in Definition 13.31 satisfy the 
following equations outside the compact set, 

S/^su = ( - |£5X|2 + ]^{Trs{CsX)fy, 

A5X' + Au~\j{CsX)i - 2u-^u'Trs{CsX) = 0, 

and asymptotic flatness requires u, X tend to 1, respectively at some decay 
rate. Using the decay conditions on u and X, we have Asu = 0(|x|~^^^^) 
and AsXi = 0{\x\~'^~'^^). As shown in |B86j . the asymptotic behavior 
implies that outside a compact set, 

u = l + -^ + Oi\x\-'-^), X' = ^ + Oi\x\-'-^), (3.3) 
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for some constants a and 6*. Note that this clearly implies that (5, vr) 
is AF-RT. Furthermore, if {M,g,TT) is AF-RT, u'"^'^{x) and (X^)"'^'^ sat- 
isfy equations with better decay: Asu"'^'^ = 0(|x|~'^"25) ^nd As{X'')"'^'^ = 
0(|x|-4-2^). Hence 

u°'^'^{x) = J^ + 0{\x\-^-^) and (X^f'^'^ix) = + 0{\x\-^-^) (3.4) 

for some vectors c, d^^^ which are quantities corresponding to the center of 
mass and angular momentum of {g,Tt). 

Theorem 3.1. [CS061 Theorem 1] Let {gij-6ij,7rij) G W'^f{M)xW^f_g{M) 
be a vacuum initial data set, where S G 1) and p > §• Given any e > 0, 
there exist ko > and a sequence of solutions {gk,Ttk) with harmonic asymp- 
totics satisfying 113. S^) . 113. 3\) so that 

h - 9k\\w-i,P(^!^j) < e, ||vr - *fcllvyif_^(A/) ^ f"^ ^ ^ ^0- 

Moreover, the mass and the linear momentum of {gk^T^k) CLfc within e of 
those of {g, vr). 

The theorem says that the solutions with harmonic asymptotics ()3.2p 
are dense among general solutions. More remarkably, the mass and the 
linear momentum which can be explicitly expressed for solutions with har- 
monic asymptotics converge to the original initial data set in these weighted 
Sobolev spaces. However, in the above theorem the center of mass does not 
seem to converge, neither is the center of mass Cj defined generally for AF 
manifolds. Therefore, we would like to modify their theorem and prove, in 
some weighted Sobolev space, solutions with harmonic asymptotics form a 
dense subset inside AF-RT solutions so that the centers of mass and the 
angular momentum converge. The precise statement is as follows: 

Theorem 3.2 (Density Theorem). Let{g-6,Tr) G W^f{M)xW^f__g{M) be 
a vacuum initial data set and [godd ^ ^even-^ ^ W^f^g{M\BR^) x w]:^^g{M\ 
Bro), where 5 G (^, 1) and p > 3. Given any e > and 5q G (0,5), there 
exist R, kQ = kQ^R), and a sequence of solutions {gk,T^k) with harmonic 
asymptotics satisfying 113. 113.^) . (f3.^| j so that {gk,Ttk) is within an e- 
neighborhood of {g, tt) in the WZf{M) X W_f_g{M) norm and 

Moreover, the mass, the linear momentum, the center of mass, and the 
angular momentum of {gk,T^k) (^fc within e of those of {g,TT). 

We first briefly describe Corvino and Schoen's construction of the approx- 
imating solutions (gfc,7ffc). 
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Sketch of the proof of Theorem fXil Let {gk,7Ck) be 2-tensors cut off from 
the original solutions {g, ir), gk = £,kg + (1 - S.k)S = 6 + ^kh, vTfc = ^fcvr where 
^fc is a smooth cut-off function 

{1 when |x| < k, 

between and 1 when k <\x\ < 2k, 
when |x| > 2k. 

S,k is chosen so that \D^k\ < f and \D'^^k\ < -§z for some constant c inde- 
pendent of k. Let 

gk = Uk9k, 

To simplify notation, we denote C-g^, by £, and we also drop the subindex 
k when it is clear from context. In order to find u near 1 and X near 
at infinity so that {g, vf) is a vacuum initial data set, we need to solve the 
following systems for u and X from the constraint equations (j3.ip 

li = u-^(^- 8AgU + (^R{g) -\tt + liX\\ + ^ {Trg{j^ + llX)f^^ = 
(divg(7f))^ = (^(div^(7r + + Au-^g^%{TT + CX)ik 

-2u'^u^Tr-g{T: + CX)^ i = 1,2,2,. (3.5) 

Consider the map r : (1, 0)+VF^'|'(M) xiy^f (M) ^ VFlf_5(M)xVF°'2%(M) 
defined by r(u, X) = (/i, divg(7f)) . It is known that Z)T(i o) is a Fredholm 
operator of index 0. For the Fredholm operator of index 0, the operator is 
injective if and only it is surjective. However, it is not clear whether Z'T(x g) 
is surjective. Corvino and Schoen enlarge the domain and utilize initial data 
sets of the form {u^'g + h, n^(7r + CX) + g) with 

$ {u^g + h, u^TT + CX) +q) = (0, 0), 

where h and q are symmetric (0, 2)-tensors with compact supports. Then 
they prove that the operator L'<I>(^ ,^) maps surjectively onto Wj2_j{M) x 
VF°f_5(M) for p > 1 and (5 e (0, 1). 

Since DT(i g) is Fredholm of index 0, we have 

W^fiM) X W^'^iM) = Ker (Z)r(i,o)) I^i, 

I^lf_5(M) X W^'^_s{M) = Range (DT^i^o)) span{yi, . . . ,Vn} 

where Wi is an A^-dimensional linear subspace and {Vi, ■ ■ ■ , Vn} is a basis 
for the cokernel of DT(^i Qy Because Z)<I>(^ ^i-) is surjective, we can choose 
{ihi,qi), ■ ■ ■ ,{hN,qN)} so that D<^(^g.^-^{hi,qi) = Vi. Supports of those 
{hi,qi) may not be compact, but there exist {hi,qi) with compact supports 
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close enough to {hi^qi) so that Vi = D^(^g.f^-j{hi, qi) still span a complement- 
ing subspace for Range (L>T(i q))- 

Let W2 = span{(/ii, gi), . . . , {h^, qw)}- W = Wi x W2 is a Banach space 
inside W^f{M) x w!'J(M) x W^f{M) x W^f_g{M). Define the map T 
from ((1,0), (0,0)) +W1XW2 to VFi'2%(M) x W^'^^^iM) by 

T((u, X), {h, q)) = <^{u^g + h, u\Tr + £X) + q). 

-Dr((i o),(o,o)) is an isomorphism by construction. Hence the inverse func- 
tion theorem asserts that T is an isomorphism from a fixed (independent 
of k) neighborhood of ((1,0), (0,0)) to a fixed neighborhood of $(5, vr). Be- 
cause (0, 0) is contained in the image when k large, there exists a unique 
{{u, X), {h, q)) within that fixed neighborhood of ((1, 0), (0, 0)) such that 
$ {u'^g + h, n2(7r + CX) + q) = (0, 0) for k large. □ 

Notice that the supports of h, q in the proof may not be uniformly bounded 
in k, but it is important in the proof of Theorem 13.21 that h, q have compact 
supports uniformly bounded in k. Therefore, we need to carefully choose 
the cokernel of LT(i g) for k large. This choice is described by the following 
lemma. 

Lemma 3.4. V, W are Banach spaces. Assume Sk ■ V ^ W is a sequence 
of Fredholm operators and Sk converges (in the operator norm) to some 
Fredhholm operator S' . IfW = RangeS" © W for some finitely dimensional 
closed subspace W , we can choose a subspace Wk in W' such that W = 
RangeS'fc ® Wk for k large. 

Proof. Since S' is Fredholm, there exists V such that V = KerS' © V. 
Consider 

by Ts{v,w) = Sv + w. Then ts,, is an isomorphism for k large since ts' is an 
isomorphism by construction and isomorphism is an open condition in the 
space of linear operators. We then have 

W = Sk{V') © W. 

Therefore, for any v £ V, Sk{v) can be decomposed uniquely into Sk{v') + 
Skiv) — Sk{v') for some v' G V and Sk{v) — Sk{v') G W. Hence v can be 
decomposed uniquely into v = v' + {v — v') as well, where Sk{v — v') G W . 
Let Uk = {u £ V : Sk{u) £ W'}. It is easy to see that C/^ is a closed space 
in V and 

v = v'e Uk. 

Note that KerSk is a finite dimensional subspace in C/^, so we can write 

V = V' ® KeiSk © Zk 
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for some closed subspace C Uk- Sk{Zk) is closed in W', so there is 
Wk C W such that Sk{Zk) eWk = W' and hence 

W = SkiV) e SkiZk) (BWk = Ranges k Wk. 

□ 

Corollary 3.5. The supports of h and q in the proof of Theorem ] 3 . 1\ can he 
chosen to he uniformly hounded for k large. 

Next, the key lemma (Lemma I3.8P used to prove Theorem 13.21 is an a 
priori type estimate for 2nd order elliptic equations Pv = f which have a 
symmetric property at infinity. Roughly speaking, if we know P and / are 
even (or odd, respectively) then we hope solutions v will be even (or odd, 
respectively) as well. However, this is not true generally because boundary 
values can affect solutions dramatically. For example, consider two harmonic 
functions ui , U2 in \ Bj^ , 

1 C ' X 

ui = - — r and U2 = — rr- 
\x\ \x\'^ 

Both |ui| and \u2\ tend to zero at infinity. However, ui is even and U2 is 
odd. They are solutions for different boundary values on the inner boundary. 
Nevertheless, in the case that the boundary value is very small, we will show 
the symmetry of the solutions is not affected much in the region away from 
the boundary. Before we state the lemma, we will give some definitions of 
the operators we consider. This class of operators is discussed in detail in 
[B86] . 

Definition 3.6. P defined as Pu = a^^{x)dfjU + V{x)diU + c{x)u is an 
elliptic operator. 



[1) P is said to he asymptotic (at rate t) to an elliptic operator P 



Pu = a'^^{x)df,u + b'^{x)diU + c{x)u, if there exist q G (3, oo), r > 0, 



and a constant C such that over the region M \ -B/j^, 

- a'^Wwl:^ + 11^' ~'^*llvF!f_, + ll'^ - c||^o,5/2 < C. (3.6) 
(2) P°'^'^ is the odd part of the operator, defined on Af \ Br^ hy 

Remark. From the definition (13.61). we have 



as i? — )• oo. 
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Definition 3.7. We say a sequence of elliptic operators P^, PkU = (^^('f:){x)9fjU+ 

Vf^^-^{x)diU + C(fc)(x)n, is asymptotic to P uniformly if given e, there exist R 
and ko such that 

+ \hk) - ^\\w!_'^i:um\b!,) ^ ^ f"'' ^ > ^0- 

Lemma 3.8. Let {Pk} be a sequence of 2nd order elliptic operators asymp- 
totic to As, the Euclidean Laplacian mM^, uniformly, and for a non-integer 
positive number a, 

Pk : VF^f (M) ^ W'Lt^liM). 

Assume {vk} C 1¥^1'^(M), {fk} C W^'^^{M) are sequences of functions. 
We also assume vf'^ G W!:f_^{M\Bno), hf'^Wwl'^iMXBn^) ^0 ask ^oo, 
and v'^ satisfy Pkv'^ = fk- Then there exist R large and ko large such 
that 

\\vf'^\\wlf^JM\2R) < ^\\fk\\wt-lt{M\BR) + ^ ■^^'^ ^ > ^0, 

where c = c{s,p,a, \\Pk - As\\op),R = R{\\Pk - ^5\\op),ko = ko{R). 

Proof. Because v'^'^'^ is not a global function, we multiply it by a smooth 
function 0, 

when \x\ < R, 
(l){x) = ^ between and 1 when R < \x\ < 2R, 

1 when |x| > 2R, 

where (j) satisfies \Dcl)\ < ^,\D'^(j)\ < for some constant c independent 
of R, where the real number R will be chosen later. We have cpv'^'^ S 
VF!f_„(M) and 



By |B861 Theorem 1.7], there exists a constant ci = ci(s,p, o) such that 

IIK''llH/!f_JM) < ci||A5(</,t;f'^)||^.-2,.(^,). 

We then consider the difference between and Pk- The operator norm 
\\Pk- As\\op-m\Br defined by 

\\Pk - Ai\\^p.M\Bn = ™P {||(i^fc - As)w\\^s-2,p^^j^^ : ||ty||vKlf_JAf) = 1> 

suppij; (Z M\Bb^ 
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tends to zero as i? — )■ oo uniformly in k because Pk is asymptotic to 
uniformly. More precisely, for u as above, 



{Pk - ^s)w\\w!_~^''^(M) - { h(fe) - I } \\D^w 



x\>R 

+C2\\b\L,o.gf,,sr,,\\Dw 



+C2 1 1 c| I ^o,,/2 (^^^^^ ^ 1 1 1 1 wl-f_^ (M) 
as i? — 7- oo uniformly in k. 



Therefore, 



U< \\wl-_JM) < - A5)(0t;r)llH^-^>p(M) + ci||Pfc(0t;r)llH/l-!-(M) 

< Ci 1 1 Pfc - A5 1 1 op;M\Br 1 1 ^'"°k'^ 1 1 (M) 



Choose R such that ci||Pfc — A^Hop-^v/^^^ < ^ and move that term to the left 
hand side. Then 

IIK'lliylf.jM) < 2ci||P,(K'')IIh/!-!-(m) 



< '^'^\\fk\\wt-'^-^{M\BR) + ^3||■yfc°''^||^y^-l.P(^^)• 



Because ||''^fc'^'^|lvi/=-i'J'(M\Bji ) — ?• as — )• 00, for e = there exists ko 
such that for all k > k^ 



alt, we have the Schauder-type estimate 

vf'^\\wlf_JM\2R) < C2||/fe|lvi/=-2.P(M\i3^) + C3 for ah k > ko 



□ 



Proof of Theorem \3.2l 

1. Estimates on and {X')"'^'^. 

We construct (^fc, Tffc) as in Theorem 13. II in the form = uf.gk + /ifc, vffc = 
u1{fck + Cg^Xk) + Recall that is the radius of which we cut off the 
original data. Again, we drop the subindex k when it is clear from context. 
By Theorem 13.11 u and X exist and satisfy the system of the constraint 
equations (13. 5p . From the constraint equation (13. 5p for u in M \ Br^^ , we 
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have 



Is ^ 2 



A^^x - - - Ivr + CX\j. + - (Tr^vr + 



- 1^ + + i(rr^(7r + CX))' 

= PlU. 

On M \ Bji^ , satisfies the fohowing equation, 

1 / 1 X""'"' 

= -i^R{g)-\^ + CX\l + -{Trgi^T + CX)j ui-x) 

- 1 / Q 



1 ( 



/1 + /2, 



where 



fi = -[-\n + CX\l + -Trg{7T + CX)j ui-x), = P^u'""' - h. 

fi contains the terms involving X. We will use a bootstrap argument to 
improve its decay rate. /2 contains the terms which have expected good 
decay rate already, such as g°'^'^ ^q^'^'"^'^ , A direct calculation tells us 

" " W]:i_^, {M\Ba^ ) ^ 1 1 /2 1 1 w'j_,{M\Bn, ) - 

We emphasize that through out this proof, c is a constant independent of k. 
From the constraint equations (|3.5p for X, we have 



(div^(£X))^ + div^(7r), + Au-^g'^^ujiT: + CX).^ 
- 2u~^UiTrg{TT + CX) = 0. 

If we compute the first term in local coordinates, we have 

(div^(£X)) . = (AgX)^ + RicikX'' 
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We define P2Xi = g^^ gSdx^ ^^"^ then P2Xi = Fi, where Fi contains the 
remainder terms from above identities, 

-divg(7r)i - Au~^g''^Uj{TT - CX).^^ + 2u'^UiTrg{T: + CX). 
Then we have 



P2 {Xi)'"''' = {FiY'''' - 

= /i,3 + /i,4 



where 



/i,3 = - (diVg(7r)i)"'''' , /j,4 = P2 {XiY'^'^ - fi^S: 



where /j^4 contains u""^"^ and X^'^'^ and we wih bootstrap to improve its decay 
rate. A straightforward calculation gives us 

\\fi'^\\w'l'^_,{M\Ba„) ^ ^' \\fi'^\\w°'^_,,(M\Bn,) ^ 
From the above, we derive the system 

Piu"'"' = h + f2, 

P2iXir'''' = n,3 + i = 1,2,3. 

We are at the stage that we can apply Lemma 13.81 (with a = 5) for each 
equation because Pi and P2 are obviously asymptotic to A5 and 

\\-^i'^'^\\w^f{M\BRg) - "^W-^iW W^f{M) ~^ as A: 00. 
Hence, there exist Ri and ki such that for all k > ki, 

W odd\\ ^ \\(-y''^\odd\\ ^ 

Once we derive these estimates, the decay rates for /i, /j^4 are improved. The 
bootstrap argument allows us to conclude that for some R2 ^ Ri, k2 > ki, 

Therefore, for any given e and 5o £ (0,5), there exist R and ko so that for 
all k > kQ, 

II odd 1 1 ^ TjSoS ^ \\( •y'''-\odd\\ ^ 

Furthermore, by Corollary l3.51 the supports of (/ifc, q^) are uniformly bounded 
in k. Hence we have 

1 1 tP^^ II c 

Similarly, we have the estimate for f . 

2. Convergence of the Center of Mass and the Angular Momentum. 
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To prove the center of mass and the angular momentum of {g, tt) converge 
to those of {g, tt), the same idea of proving convergence of the mass and the 
Hnear momentum in [CS06j is employed. 



\Cf{g)-Cn9)\< 



(9)- [ {Ri^-\R,g^,){Y'^Yyida, 

J\x\=r ^ 



+ 



+ 



J\x\=r ^ J\x\=r ^ 



C?{-g) 



The first and the third terms can be written as integrals over > r} by 
the divergence theorem. We will also use the fact that {|x| > r} is centrically 
symmetric, i.e. x, —x G {|x| > r}, to estimate those integrals over {|x| > r}. 
For the first integral, 

Cf{g)- [ 

J\x\=r ^ 

= [ x''Rg-2x''{g'^ -5'^){Ri,-\Rgg,j) 

J\x" 



'|x|>r 



+ {R,,--R,g,,){Y'^YTU^d^o\g. 



(3.7) 



Let H = {x ^ M \ Br : > 0} he the half space. Then we can rewrite the 
above integral as follows: 



/ x^{Rf''^g + Rg{-x)^g'""')dx 
Jh 

2x'^{g'^r''\R,j-\Rgg,,)^dx 



l^2x'^{g'^ - 5'^) - \R°fg,, - }-R^{-x)g^'')^g + {R,, - 



dx 



+ j^{Y'^Y(R:if - ^Rf% - ^-Rgg'^^'^ytig'^ dx 



We substitute for Rg in the first integral by using the constraint equation 
Rg = — ^(Tr^vr)^ + |7r|^. We then bound above integrals symbolically by 
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the following, 

ci / \x\(\Tr^''^'^\\n\ + \TT\'^\g°'^'^\ + \g°'^'^\\D^g\ + \g-5\\D^{g°^^^^^ 
+C1 ^ \x\^(\D\g"'"')\\Dg\ + \D^ g\\D{g"'''')\ + \D^g\\DgU'"'\) dx 

^ ^ I ll^e'^e"-!! Il^ll i ll^l|2 ||„odd|| 



MW''\\w-^. Mw-M-^\ (3.8) 

—1—0 —o / 

where ci,C2 are constants independent of g and vr. The weighted Sobolev 
norms above are over the region > r}. Similarly, 

Cf{g)- [ {R,,-lR-g,,){Y'^)\lda-, 

J\x\=r 



^ „ / 1 1 - even \\ iI|;=l||2 ii -odd \ \ 



-1-5 



+ Mw^-y-'- (3.9) 



For the surface integral, we can assume k ^ r (recall k is the radius of 
which we cut off the original data) and then g = u^g on {|x| < 2r}. Hence 
we have 



J\x\=r ^ J\x\=r ^ 

= [ ({Rij - u^Rij) -\{Rg- u'^mgij) (y")V^' dag 

J\x\=r ^ Z / 

= {a-u^)RiJ-la-n^)Rg9^J)iYn'l^i,dag 

+ / (u\Rij - R^j) - lu^iRg - R)giA{Y''yu^g dag 

J\x\=r ^ Z / 

< C6 max (|1 - max (\D'^g\\x\'^yr'^-'^^ 

\x\=r V / |x|=r V / 

+C6 max \D^{g - 5)] jxp+^r^"''. 

\x\=r 

Recall the Sobolev inequality for weighted Sobolev spaces (see [B86J). For 
n — kp < where n = dimM = 3 , we have ||u)||^o,oo^j^^^^ < c||u)||p^fc,pj.^^^^. 

Therefore, for p > 3, 

max\D^{g -g)\\x\^+^ < maic \D^{g - g)\\x\^+^ 

\x\=r M\Br 

= \\D^{9 - Mw'^jZ.iMXB,.) ^ c||I?'(5 - 9)\\wy,_,(M\B,.) 

< 49 - 9\\w^f^M\B,y 
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Similarly, using the Sobolev inequality on the other term, we derive 

J\x\=r ^ J\x\=r ^ 

To conclude the argument, for given e, let Ri be a constant so that in (j3.8p . 
max{c2(||7r--||^i,p_J|7r||^i,._^ + ||7rf^,,,_J|g°^^||^2,p_^ + ||/'^''||^2,._J|5||^ 

For this e, there exist R2 > -Ri and /cq so that 

in II ^ e II - II ^ ^ 



„ - 1 1 ^ 1 1 -odd II ^ 



•-2 



As a result, using the estimates (|3.8|) (j3.9|) (|3.1Up and try r in those estimates 
equal to R2, we can conclude for all k > ko, 

\Cfig)-Cf{-gk)\<e. 

To prove that angular momentum of g is close to that of 5, we notice that 

Ja- [ 7rjk{Zyu^gdag= [ dWg{7Tjk{Zy) dVOlg 
J\x\=r J \x\>r 

The first term is zero since irjk is a symmetric tensor and is a Killing 
vector field. The other terms after integration are bounded by Cr^^^^. Then 
the rest of the argument works the same as the case of center of mass. □ 

In the proof of convergence of the center of mass and the angular momen- 
tum, we showed that the limits at infinity can be approximated by surface 
integrals at the finite radius and the difference for surface integrals at the fi- 
nite radius is arbitrarily small. Since we will use this argument several times 
through this article, we formulate the argument into the following lemma. 

Lemma 3.9. Let F{g) he a vector field depending on x , gij , D gij , D"^ gij 
smoothly defined in M \ Brq. Let a,b > 0, and let r be the radius. As- 
sume the following estimates hold for k > r, and for g,gk satisfying the 
assumptions in Theorem \3.2\ (Density Theorem). 
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(1) 



J\x\>r 



d.\YgF{g) dvolc 



< c k' 



+ m 



odd I 



W 



if_,(M\B.)"^ llw!f_,(M\S.) 



(2) 



I 1 1 odd II n il I II „||2 

^ \\w^f_s{M\Br)\\9\\w^f{M\Br) \\y\\w!_f(M\Br) 



/ -f'Cf ) • t^g d(Tg - / F{gk) ■ Vg^ dug^ 

J\x\=r J\x\=r 



T-\ (3.11) 



< c( 111 - '"l|w^3,^,^^^^g^) 



+ Il5 - 9k\\w^^v_^iM\Br) V^- (3-12) 



Then given e > 0, there exists ko such that 

/ F{g) ■ Ug dag = lim / F{gk) ■ Ug dag + e, for all k > ko. 

J\x\=r ^^°°J\x\=r 



lim 



Remark. If we replace the normal vectors and the volume forms of the 
integrals in the above assumptions by those with respect to the induced metric 
in the Euclidean space, the analogous result is the following: assume the 
following estimates hold 



(1) 



J\x\>r 



divsF{g) dvolo 



Wt'i_,iM\Br) II " \\wif_,iM\Br) 



+ m 



odd I 



W^f_,{M\Bry^9 \\w^'f_giM\Br) 



+ \\9"'^'^\\w^f_^iM\Br)y\\w^_f(M\Br) + y\\w^'P{M\Br))^ (^•^^) 



(2) 



F{g) ■ vq dao - / F{gk) ■ uq dao 

\x\=r J\x\=r 



<c(||l «l|vi/3'P(M\B^)ll5 5'fc|lvi/3,P_^(Af\B^) 1^^- 



(3.14) 



Then given e > 0, there exists ko such that 
lim / F(g) ■ vq dao = li™ / P{9k) ■ dao + for all k > kg. 
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4. Properties of the Intrinsics Definition 



In [B86| Section 4], Bartnik proves that the ADM mass is a geometric 
invariant and satisfies natural properties. Along the same lines , we will 
show that the intrinsic definition of center of mass ()1.4p is well-defined and 
has corresponding change of coordinate under the transformation at the 
infinity. We will first show that the intrinsic definition is robust in the sense 
that we can integrate over a general class of surfaces, and those integrals 
converge to the same vector. 

Proposition 4.1. Suppose {M,g,K) is AF-RT Let {Z^fcl^i C M 

he closed sets such that the sets Sk = dD^ are connected 2 -dimensional 
-submanifolds without boundary which satisfy 

rk = inf{|x| : x G Sk} — )• oo as k ^ oo, (4-1) 

r^^area(S'fc) is bounded as k ^ oo, (4-2) 

vol{Dk\D-} = 0{rl), 

where = Dk D {—Dk} and 3~ is a number less than 3 . 

(4.3) 

Then the center of mass defined by 
is independent of the sequence {S^}- 

Remark. The first two conditions (|4.ip . (|4.2p on Sk are the conditions con- 
sidered by Bartnik in [B86j to ensure the ADM mass is well-defined. The 
volume growth condition (|4.3p allows us to consider a general class of surfaces 
which are roughly symmetric; that is, the non-symmetric region Dk \ of 
Dk has the volume growth slightly less than the volume growth of arbitrary 
regions in M. 



Proof. As the equation (|3.7p . the divergence theorem gives us 

[Rij - \Rg9i3)Y{a)^l dag = {Rij - ]^Rg9ij)Yl^-^vl dag 

We can decompose the integral over Dk \ Di = {Dk \ D^} U {D^ \ D^} U 
{D^ \ Di} into three integrals, where D^ = Di U {~Di}. The first integral 
over {Dk\D^} converges because of the volume growth condition ()4.3p . The 
second integral over {D^ \ Df} converges because the region is centrically 
symmetric and the initial data set is AF-RT. After taking limits, the right 
hand side has a limit independent of the the sequence {Sk}- 
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□ 

Assume {x} and {y} are two AF coordinates on (M \ g). Assume F 
is the transition function between these two coordinates and y = F{x). In 
|B86l Corohary 3.2], it is shown that the only possible coordinates changes 
for AF manifolds at infinity are rotation and translation. More precisely, 
there is a rigid motion of M^ (Oj, a) G 0(3, M) x M so that 

\F{x) -{Ox + a)\e W^'°°{M.^ \Bii,). 

Similarly, if we use the same argument in that corollary for AF-RT mani- 
folds, we derive 

F°'^'^{x) G W^'^iR^\BR,). 

Because the center of mass is a quantity depending on the coordinates, 
we now show that centers of mass in {x} and {y} coordinates have the 
corresponding translation and rotation. An interesting phenomenon is that 
compared to rotation, translation is a more subtle rigid motion. If the 
translation a is not zero. Density Theorem is involved in the proof. 

Theorem 4.1. Let {x} and {y} be two distinct AF-RT U.3\} coordinates for 
{M,g, K) satisfying the change of coordinates as we described above. Assume 
Ci^x in-d Ci^y are the centers of mass defined by the intrinsic definition 1^1-4^ 
in these two coordinates, then 

Ci,y = OCj^x + a- 

Proof. The metric g in the cotangent spaces induced from {x} and {y} can 
be written locally as 

ds^ = gki{x)dx''dx' = gijiy)dy'dy^ = gij{F{x))d{Olx'' + a)d{Ojx' + a) + 0{\x\-^] 
Expanding the above terms, we get 

gki{x) = gij{F{x))OlOi + ei{x) 

where the error term ei G M^^'^ and ef^'^ G W]_'^. Then a straightforward 
calculation gives us 

and the formulas for Ricci and scalar curvatures, 

Rki{x) = Rij{F{x))OiO'i + esix) 

R{x) = R{F{x)) +e4(x) 

where eq{x) G 0(|x|~^) and e°'^'^(x) = 0(|x|~^) for q = 2,3,4. To calculate 
Ci^y, Proposition 14.11 allows us to integrate over {\y — a\ = r}, 

^ty = TTT— 1™ / (R^Jiy)-lRiy)g^J{y)) (|yp<5°^-2y-y') dao. 

levrm r^oo V 2 / V ^ \y - a\ 
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We would like replace y = F{x) in the above identity. First we have 

|y|2^°i - 2y"y' = {\xf - 201x1^01x^)6''' + 20x ■ a6"' 
-2{a''Oix^ + O'^x'^a') + 0(1), 

T 1 = + 

\y ~ Q^l fI 

where eg G W^'^ and ef'^ G VF^'g"". 

J|j;-a|=r ^ ^ ^ ^ ^ |y - a| 

J|x|=r ^ ^ ' ^ l-^l ' 

+ / - \R{F{x))g,,{F{x))) hOx ■ aS'^' - 20^x^a*)) ^ ddo 

J\x\=r ^ ^ ' ^ Pl 

+ 0(r-i). 

We claim that the first integral 1\ converges to IGvrmO^Cj ^, the second 
integral I2 converges to IGvrma", and the third integral I-^ converges to 0. 

h = I (Rij{F{x)) - \R{F{x))g,j{F{x)))oioU\x\'0'^-20''^x^x'^)^dao 

J\x\=r ^ ^ ^ ^ ^ l-^l 

= ^ ^_ ^Rkiix) - ^Rg{x)gki{x)) (1x^6^" - 2x''x'=) ^ dao + 0{r~') 
= 167rmO'1Cl^ + 0{r-^). 
By using the formula for ADM mass (jl.Sp . 

f 1 ' 

I2 = a" ^ ^ (i?H(x) - -Rg{x)gki{x)') ( - 2a;'=^) d(Jo + O(r-i) 

= levrma" + 0(r~^). 

To see that Is converges to 0, we simplify the expression by letting b = O^a, 
and then 



h = o 



\x\=r 



Rki{x) - \Rg{x)gui{x)) {2x ■ h5^^ - 2x^6'=) ^ da^. 
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By the fact that (2a; • b5^^ — 2x^6^) ■ ^ = 0, the scalar curvature term 
vanishes after taking a hmit, so we only need to prove 

J \x\=r ^ ^ l-^l 

It is sufficient to prove 

[ Rki{x)(2x ■ h5^^ - Ix^b^^vldug = O(r-i). 

J \x\=r 

A straightforward calculation shows this identity holds if the metric is con- 
formally flat outside a compact set, so we would like to approximate g by 
g which have harmonic asymptotics, and then apply Lemma l3.9i We need 
to check if conditions (j3.1ip and ()3.12p hold. Using the divergence theorem, 
we have 

diYg(Rki{x){2x ■ b5^^ - 2x^h^)) dvol^ 
{diYgRki){2x ■ h6^^ - 2x^h^) + Rki{2l}5^^ - 25^^h^) dwolg 

\x\>r ^ 



sum to =0 



+ / Rki [2x ■ bS^'' - 2xH^)g'^TZn ^volg. 

J\x\>r 

Symbolically, the above integral is bounded by 

divg(^i?fc;(2;)(2x • U^^^ - 2x^b^)^ dvolg 



^^ni II \\^even\\ i IUI|2 

Hence condition (jS.lip holds. The condition ()3.12p holds because the Sobolev 
inequality implies 

/ Rki{x) hx ■ b5^^ - 2x^b^) v\ dag - j Rki{x) hx ■ b5^^ - 2x^b'') v 

J\x\=r J\x\=T 



By Lemma 13.91 we conclude 

lim I Ru{x)i2x-bb^^ -2x^b^\v\dag 

r-'°°J\x\=r ^ ) y ^ 

= ^lim y" Ru (x) (2x ■ bb^^ - 2xH^^ v\ da-g + e = e. 



Since e can be chosen arbitrarily small, converges to 0. Therefore, 

1 



Cty = Am h+h+h = O^Cl + a" 



□ 
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5. The Corvino-Schoen Center of Mass 

In ]CS06] . Corvino and Schoen define the center of mass (jl.Sp for AF- 
RT manifolds. In this section, we will show that the intrinsic definition we 
propose is actually equal to the Corvino-Schoen definition (Theorem [1]). In 
other words, the intrinsic definition (|1.4p is a coordinate- free expression of 
the Corvino-Schoen definition. 

Assume g is the approximating solution in Theorem 13.21 (Density Theo- 
rem). A straightforward calculation shows Ci{g) = Ccsid)- As proven in 
Theorem l3.2t Cj{g) can be approximated by Cj{g). If we prove that Ccs{g) 
can also be approximated by Ccsig)-, then Theorem [T] follows. 



Proof of Theorem 1. First notice that the Corvino-Schoen definition is equal 
to the following expression where the normal vector and the volume form 
are with respect to the induced metric in the Euclidean space 



n 



1 



CS-T7, I / x'^iOij,!- 911,3)^1 dao 



A direct calculation shows 



|a;|=r 



x"igij,i - gii,j)i^odao - / [hia^l - hiiu^) dao 

x\=T2 J \x\=r2 



x\=ri 



{hiai^' - HuUq) dao 



x\=ri 



ri<\x\<r2 



x°'{gij,ij - gujj) dvolo 



It is easy to see that the conditions (jS.lSp and (j3.14p in Lemma 13.91 (and 
the remark afterward) hold because gij,ij — gujj is the leading order term of 
Rg, and by the constraint equation ()3.ip . 



x°'i9ij,ij - gujj) dvolo 



x\>r 



.odd\\ ^1-25 



and 



x°'{gij,i - gii,j)^o d<^o - / {hiavl - hiiVQ) dao 



x°'igij,i - gu,j)i^odao - / {hiavl - hiiVQ) dao 

\x\=r J\x\=r 
2-<5 
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Then we can apply Lemma 13.91 and have 



^^°°J\x\=r J\x\=r 

= hm / x"{gij^i - gii,j)vi dao - / {hiai'l - hii^o) dao + e 

r^'^J\x\=r J\x\=r 



'\x\=r J\x\=r 



□ 



6. The Huisken-Yau Center of Mass 

In the case that (M, g, K) is SAF (US]), Huisken-Yau |HY96) and Ye |Y96] 
prove the existence and uniqueness of the constant mean curvature fohation 
in the exterior region of M, if m > 0. Huisken and Yau use the volume 
preserving mean curvature flow to evolve each Euclidean sphere centered 
at the origin with a large radius and they show that the Euclidean sphere 
converges to a surface with constant mean curvature. Furthermore, they 
prove those surfaces with constant mean curvature are approximately round, 
and their approximate centers converge as follows: 

r, r dcTn 

C^y = lim ^^^^^^ — , a = 1,2, 3 (1.7) 

where {Mr} are leaves of the foliation and z is the position vector. They 
define this to be the center of mass for M. 

Instead of using volume preserving mean curvature flow, Ye uses the im- 
plicit function theorem to find a surface with constant mean curvature which 
is perturbed in the normal direction away from the Euclidean sphere Sji{p) 
for some suitable center p and for the radius R large. Although it is not 
emphasized in his paper, p can represent the center for each constant mean 
curvature surface because this surface is roughly round. Furthermore, we 
will show that p converges to Ccs and use this fact to prove Chy = Ccs- 
Then C/ = Chy follows because the Corvino-Schoen center of mass Ccs is 
equal to the intrinsic definition Cj by Theorem [TJ Before we prove Theo- 
rem [21 we will need a technical lemma (Lemma 16. ip which suggests that the 
center of the surface p converges to Ccs- 

Let z/ be the outward unit normal vector field on 5/j(p) with respect to 
the metric g. We parallel transport along the normal direction so that it is 
well-defined in a tubular neighborhood of Sii{p). Hence the mean curvature 
on y G Sji{p) is div5^(p)Z/ = div gU, the divergence operator of the ambient 
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manifold M. 

^\y-p\ 

m 3m? 1 {y'^ — p'^){y^ — p^)\y^ — d 

A straightforward calculation gives us the mean curvature on Sji{p) is equal 
to 

^. 2 4m 6m{y — p) • p 9m? 

^ \y-p\ |y-pP \y-p\'^ \y-p\^ 

+ 2^^^i'^^2/j |^_^|3 +^P^J[y) |^_^|3 

where |^o| < iF(l + |p|) for some constant c depending only on the metric 
9- 



Lemma 6.1. For R large, 



-SmirC^s + OiR'^), a = 1,2, 3. (6.2) 



Proof of Lemma \6.1[ We denote the first integral by 

J\y-p\=R V2 

In the proof, we will rewrite I{R) , and then some cancellation allows us to 
rearrange the left hand side of (j6.2p so that it has an expression correspond- 
ing to egg. 
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Since the coordinate is only defined outside a compact set of M, we can use 
the divergence theorem only in the annular region A = {R < \y — p\ < Ri}, 



Pij,i . [2 '^volo 



2Ja'^'"\ \y-p\ 



,k 



+ 2 [ J ^ d^olo. 

Using the divergence theorem and simplifying the expression, we get an 
identity containing purely boundary terms 



l{Ri) = I{R) + J{Ri) - J{R), for all Ri > R 



where 



2 7|y_p|=K \y-p\ 
J\y~v\=R \y-pr 

1 r yOi _ pa ^ r — pi- 

+ 7T / Pin r"'^0 + 7^ / Pia-, rdaQ. 

2 J\y-p\=R \y-p\ 2 J\y-p\=B. \y-p\ 

To prove that I{R) = J{R), we would like to apply Lemma |3.9[ It is easy 
to check that the conditions (j3.13p and (|3.14p hold, so we get 

Z{R)-J{R)= lim Z{Ri)- J{Ri)= lim Zp{Ri) - Jp{Ri) + e, 

Ri-^oo -Ri— >-oo 

where Ip and Jp denote the integrals that we obtain by replacing pij by pij 
in X and J', where pij are 0(|y|~^)-terms in the approximating solutions gij 
in Theorem 13.21 (Density Theorem) . More precisely, gij has this expansion 



and Pij is defined by 



gij = [l + — )6ij + 0{\y\ 



Pij = gij - (l + 
Because p°j'^ = ^ + 0(|y|^^) as in (j3.4|) . it is easy to see that 

lim Ip{Ri) - Jp{Ri) = 0. 
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Therefore, we conclude X(i?) = J{R). We then replace X{R) by J{R) in 
the identity (|6.2p and derive that the left hand side is equal to 

iff -p' 

^ \J\y-p\=R \y-P\ 

Pia-, r - Pa-, r dao ■ 

y-p\=R \X-P\ \X-P\ J 

We rewrite the above integrals into the summation of an expression of center 
of mass (jl.Sp and the remainder terms. Because the explicit calculations 

2m' 



using gij — Pij = i^j ^ij Sive us the remainder terms has lower order, 
we derive 

iff y^ -p' 

~o / - P"){9ij,i - 9ii,j)i r dcro 

^ \j\v-v\=R \y ~ p\ 

y1 — pi- yO! — pOi \ 

9ia-, r - gu-, r dao 

y^p\=R \x-p\ \x-p\ J 

iff y^ -p' \ 

+ / ~ ((s-ij,* - Pij,i) - i9ii,j - Piij)) 1 r dao 

^ v\y-p\=ii \y-p\ J 



1 f .y -p , .y -p , \ 

"o / [9ia-Pia)] r - {9ii - Pii)-^ i-dao] 

2 yJiy_pi=n \x-p\ \x-p\ J 

-87rmCgs + 0{R-^). 

□ 



Proof of Theorem\M Let Fp^R : 5i(0) — )• M be an embedding defined by 
y = Fp,r{x) = Rx + p, that is, Fp^/j(5i(0)) = Sr{p), the Euclidean sphere 
centered at p with the radius in M. We consider the perturbation along 
the normal direction on Sr{p) defined by S = {y + A(/)z^ : y G Sr{p)} for 
a parameter A > 0, and for cp G C'^'"(>S'/j(p)) with ||(/)||c'2,q < 1. We denote 
the mean curvature on S by H{p, R, Xcj)). Using this notation, H{p, R, 0) = 
divgi^, the mean curvature of Sr{p). By Taylor's theorem for mappings 
between two Banach Spaces, we have the following expansion in the 0- 
component at 0, 

H{p,R,X(l)) = H{p,R,0) + dH{p,R,0)X(j) 

+ {I- s){(fH{p,R,s{\(t)))\^(t)(k)ds 

where dH and d^H are the first and second Frechet derivatives in the 
component. In our case, dH(p, R, 0) is the linearized mean curvature oper- 
ator on Sr(j)), i.e. 

dH{p, R, 0) = A5^(p) + \A\l + Ricii^, v), 
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where A^^^i^p) is the Laplacian operator on Sr{p) with respect to the induced 
metric from g, A is the second fundamental form on Sr{p) and Ric{-,-) is 
the Ricci curvature of M. In [HY96j . the estimates on the eigenvalues Ai, A2 
of A are derived and 

\A\l = Xl + Xl + 0(i?-3) = A + OiR-^), Ric{u, u) = 0{R-'). 

For the second Frechet derivative in the Taylor expansion, we have 

d^H{p, R, s{\<j)))\^<j)^ = ^H{p, R, t{\ct>)) 

The right hand side is the second derivative of the mean curvature of the 
surface {y + s{X(p)i' : y £ Sr{p)}. For R large, the unit outward normal 
vector field on {y + s(A(/>)z^ : y £ Sji(j))} is close to v, and a straightforward 
calculation gives us 

Aj^ x^\Dm x^m ^ 

R3 i?2 + ^ 



^H{p,R,t{X4>)) 



Y /V -L^r r \3| -11 n ill n2 ii 1 



where the constant c is independent of p, R, (p. 

To manage those terms in the Taylor expansion of the mean curvature for 

H{p,R,X<l)) = H{p,R,0) + As^(^p)X(l)+{\A\l + Ric{i^,i^))Xcl) 

+ {1- s)(^d^H{p,R,s{X(p))X'^^^^ds, (6.3) 

let G and Ei be defined as follows where G is the lower order terms of the 
mean curvature of Sji{p) from (16. ip . 

1 {y' -p'){y^ -p'){y'' -p^) ,„ {y' - p')iy^ - p') 
G(y) = -Pij,k , f3 h2pij ■ ^3 

2 \y — p\ \y — p\ 

_ y^ - p' Pa , 1 y^ - p' 

Pij,i I II I "T nPiij I I 1 

\y -p\ \y-p\ 2 \y -p\ 

and 

Ei{y) = Eq + {\A\l - ^^X<P + Ric{v,v)X^ + {I - s)[d'^H{p,R,sX^))X^(t)(P^ds. 

G is the good term that gives us Gcs as indicated in Lemma [6TT] and Ei{y) is 
an error term bounded as follows for some constant c independent of p, R, (p. 



\Ei\ < — (l+|p|)+— (^A|<A|+i?A^|2?</.||</.|+i?^A^|D0r+i?^A^|</.||Z)0||DV| 
From identities (|6.1|) and (j6.3|) . 

XT/ , , X 2 4m 6m(y — p) • p Om^ ^/ \ , a , 2 , , ^ , , 
H{p,R,Xct>) = -^-^+ ^ +-^+G(j/)+AA5^(rt0+j^A</>+i?i(y). 
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To find the surface with constant mean curvature, we need to find p, R, 4> so 
that 

2 Am 

H{p,R,X^) = --^. 

It is equivalent to solving 

= + ^ + Giy) + AA,^(,)0 + -|a0 + E, (y) (6.4) 

We pull back the equation ()6.4p via the map -Fp,_R, and we get the following 
equation on 51(0), 

= _^ + -^ + GoFp,jj(x) + AAs,(o)V'(a:) + ^AV'(x) 

+Ei o Fp^nix) (6.5) 

for p,R,Tp{x), where A5j(o) is the Laplacian on 5*1(0) with respect to the 
pull back metric and ip = (p o Fp^ji{x) = (f){Rx + p) is the pull back of (/>. 
Define the operator L by 

L = -Ao - 2 : C2'°(5i(0)) ^ C°'"(5i(0)) 

where Aq is Laplacian on 5'i(0) with respect to the round metric induced 
by the Euclidean metric on M^. Because the metric g is asymptotically flat, 
the difference between Aq and A5^(o) is small and can be treated as a error 
term. Therefore, the identity (|6.5p is equal to 

6mx • p Om^ ^ ^ ^ ^ 

= ^3 + + Go Fp^nix) - -^ALV'(x) + o Fp,i?,(x), 

where Ei has the same bound as Ei. We let A = for some fixed a G (0, 1) 
and multiply R'^'^"' on both sides of the above equation, then 

Furthermore, since D^^ = (Dy(j))R, 

|i?ioi^P,ij(^)l < 7^4(i + N) + ^(^;r^ + ^^2^ + ^ + 

In order to find p, R and Tp to solve (|6.6p . first we perturb p = p(i?, V') so 
that the right hand side of (16. 6p is inside RangeL for any R and ip. We will 
also show that p = Cos + e where e is the error term containing lower order 
terms in R and WipW- Second, using an iteration process and the Schauder 
estimate, we can find a solution ip for R large. 

1. Perturb the Center p. 

L has a kernel equal to spanjx^, x^, x^} because translation preserves the 
mean curvature. Since L is self-adjoint, C^'°'{Si{0)) has the L^-orthogonal 
decomposition C'^'"(5i(0)) = RangeL © spanjx^, x^, x^}. We would like to 
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find p so that the right hand side of (|6.6p is ortliogonal to span{x^, x^, x^}. 
That is, we want to find p so that for a = 1,2, 3, 

f „ fQmx ■ p ^m? \ 

+ [ o Fp^Rix) + R^+'^Ei o Fp^Rix)) dao = 0. (6.7) 

-'5i(0) ^ ^ 

We calculate each term above separately. A direct calculation gives the first 
integral 

f „ / 6mx ■ p 9ni^ \ Smirp" 
From the area formula and Lemma l6.ll we have 



x°i?2+-G o Fp^Rix) dao = [ yl-P-R^+-G{y)R-^ dao 

Si(0) J Snip) ^ 



/ {y--p-)G{y)da, = ^^^^^Oi^R 
JSr{p) ^ 

Moreover, the error term can be bounded by 
/ x''R^+^EioFp,R{x)da^ <^(l + |p|) 



Since the ADM mass m > 0, we can choose p 

p{R,ij) = Ccs + e{R,ij) (6.8) 
so that the identity ()6.7p holds, where 

\e{R,n < |(1 + \P\) + + + W? + mD^\\D^^\). 

2. Find the Solution ijj by Iteration. 

We consider the isomorphism L : (KerL)^ — )• Range(L) for (KerL)-*- C 
C72'°(5i(0)) and Range(L) C C70'°(5i(0)). If we denote the right hand side 
oi^hy f{p,R,il^),\.e. 



2 



f{p, R, ^) = - 1^ - ^'^"G- o Fp^R{x) - R^+-E, o FpM^), 

we know f[p{R, ip), R, ip) G Range(L) for any R, ip with ||^/;||(^2,ti < 1. There- 
fore, any ^|)Q with ||^/>o||(72." ^ li there is V'l £ (KerL)"*- such that 

Li,,=f{p{R,^|;o),R,^|^o)■ 
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Moreover, we use the Schauder estimate and the fact that ■01 £ (KerL)-*", 
||V'i|lc2.- < c||/(p(ii,V'o),^,V'o)||c-o,« 

^ -^a^\ccs\ + 1) + -^sm\c^.'^ + m\W) 

For any R large enough (independent of ipo), we have HV'i Hc^.q ^ 1- We 
continue the iteration process and get a sequence of functions {V'fclfcLo 
isfying 

Lipk+i = f{p{R,ipk),R,'4'k) and \\^pk+l\\c^,<^ < I- 

The Arzela-Ascoh theorem says that there exists t/joo G C*^'^ (5*1(0)) such that 
a subsequence of {ipk} converges to -000 in C^''^ for < /i < a. Moreover 
■000 is a solution to (|6.6p . 



^■000 = f{piR,tpoo),R,1poo)- 

Let (/>oo(y) = -000 o Fp^liiy), then the surface M/j = {z : z = y + ^i/, y G 
S'ij(p)} has constant mean curvature equal to ^ — 

To complete the proof of Theorem [21 we need to compute 



lim ^ ■ (6.9) 

By the uniqueness of the constant mean curvature foliation, the {Mr} are 
equal to those constructed in |HY96j . and therefore (I6.9p converges to the 
Huisken-Yau center of mass C^y- We now prove that (|6.9|) also converges 
to the Corvino-Schoen center of mass C^^. Let F be the diffeomorphism 
defined by F{y) = y + ^z^, then 

where JF is the Jacobian from the area formula, JF = 1 + 0(i?~^~"). Now 
we can use the fact that the area of Euclidean sphere is 0{R?) and the 
estimate for the center p in ()6.8p to conclude 

= C^s + <R^<t^oo) + 0{R-''). 

Therefore, after taking limits, the Huisken-Yau center of mass C^y is equal 
to the Corvino-Schoen center of mass C^g and, therefore, is equal to Cf. □ 
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ON THE CENTER OF MASS OF ISOLATED SYSTEMS 
WITH GENERAL ASYMPTOTICS 



LAN-HSUAN HUANG 



Abstract. We propose a definition of center of mass for asymptot- 
ically flat manifolds satisfying Regge-Teitelboim condition at inflnity. 
This definition has a coordinate- free expression and natural properties. 
Furthermore, we prove that our definition is consistent both with the 
one proposed by Corvino and Schoen and another by Huisken and Yau. 
The main tool is a new density theorem for data satisfying the Regge- 
Teitelboim condition. 



1. Introduction 



A 3-manifold M with a Riemannian metric g and a two-tensor K is called 
a vacuum initial data set {M,g,K) if g and K satisfy the constraint equa- 
tions 

Rg-\K\l + H^ = 0, 
dwg{K) -dH = 0, (1.1) 

where Rg is the scalar curvature of M and H = Trg{K) = g^^Kij. We say 
(M, g, K) is asymptotically flat (AF) if it is a vacuum initial data set and 
there exists a coordinate {x} outside a compact set, say B^^, in M such 
that, for some 5 G (1/2, 1], 

gij{x) = 6ij + hij{x), hij = 0{\x\"^) Kij{x) = 0{\x\-^-^) 

9ij,k{x) = 0{\x\-^-^) Ki^^kix) = 0(1x1-2-^) 

9iM^) = 0{\x\-^-^) Kij,ki{x) = 0(1x1-3-^) (1.2) 

and similarly for higher derivatives, up to the third derivatives on g and up 
to the second derivatives on K. For AF manifolds, the ADM mass m is 
defined by 

m = lim / iqij i — qa j) ul dan-, 

where {|x| = r} is the Euclidean sphere, Vg is the unit outward normal vector 
field with respect to the metric g, and dag is the volume form induced from 
(M, g). We remark that if one replaces Vg by the unit normal vector field vq 
with respect to the Euclidean metric and dag by the volume form dao with 
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respect to the induced metric from the Euchdean space, the hmit is in the 
same. Similarly, the remark holds for ()1.4p . (jl.Sp and (jl.Sp . 

For the case in which we are interested, we require asymptotic symmetry 
on (M, g, K). We say that (M, g, K) is asymptoticahy flat satisfying Regge- 
Teitelboim condition (AF-RT) if (M, g, K) is AF and g, K satisfy these 
asymptotically even/odd conditions 

g°4\x) = 0{\x\~^-') K^p^ix) = 0{\x\-'-') 
{g?f')^,ix) = Oi\x\-'-') {K^r),ki^) = Oi\x\-'-') (1.3) 

and on higher derivatives, where f°'^'^{x) = f{x) — f{—x) and f™^^{x) = 
fix) + fi-x), |RT74j . Notice that f"'^'^ and /^^^" are only defined outside 
Bfig in which the coordinates {x} are defined. It is proved by Corvino and 
Schoen in [CS06j that AF-RT manifolds form a dense subset of AF manifolds 
in some suitable weighted Sobolev spaces. 

For (M, 0, K) satisfying AF-RT, we propose an intrinsic definition of cen- 
ter of masqj suggested by Richard Schoen, 

= iLl - ^^^^-) ^^)^^^"- « = ^' 2' ^ ^'-'^ 

where Rij is the Ricci curvature of M and Y(^a) = (|2;p5"* — 2x"x*)^ is 
a Euclidean conformal Killing vector field. We use the Einstein summation 
convention and sum over repeated indices; though sometimes we employ 
summation symbols for clarity. The definition is intrinsic because it is a flux 
integral at infinity of the three-dimensional Einstein tensor contracted with 
two vector fields y(a) and Ug, and Y^^) has a geometric meaning. Moreover, 
the surface of integration can be defined more geometrically (proposition 
13. ip . If {M,g,K) is AF-RT, the limit converges. Note that the above ex- 
pression is not defined when the ADM mass m is zero. However, the center 
of mass can not be well-defined when m is zero because the positive mass 
theorem says that M is actually the Euclidean space, if the scalar curvature 
Rg>0 [SY791ISY8T] . 

This intrinsic definition is motivated by a similar expression of the ADM 
mass when Y^^k) in definition (II. 4p is replaced by {— 2x*^}, the radial di- 
rection Euclidean conformal Killing vector field, 

m = lim / (Rij-l-Rggijj{-2x')u^gdag. (1.5) 



The intrinsic definition can be generalized to n-dimensional manifolds if we replace 
the factor by a suitable constant depending on n and ujn, where aj„ is the volume of 

the unit ball in R". Hence the following arguments work more generally for n-dimensional 
manifolds, but for simplicity, we assume n = 3. 
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The Euclidean conformal Killing vector fields {— 2x*^} and Y(^) generate 
dilation and translation near infinity. A detailed discussion about these two 
vector fields can be found in |C W08j . Another motivation comes from the 
spatial Schwarzschild metric 



9''ix)=[l + 1 6 = 1 + — + —-j^ + + 0(|x|-3^ 

\ 2\x — p\ J \x\ \x\'^ 2\x[^ 

If we replace the metric in (jl.4p by , Cj is precisely the vector p which 
indicates the center of the manifold. It is worth mentioning that although 
in the Schwarzschild case p is not a point in the manifold, one can develop 
polar coordinates using concentric spheres centered at p. 

Other notions of center of mass have been proposed. Huisken and Yau 
[HY96| defined the center of mass for (M, g, K) which are spherically asymp- 
totically flat (SAF), i.e. {M,g,K) is AF and 

9ij{x) = (^^ + +Pij (1-6) 

Pijix) = 0(|x|-2), = 0(1x1-2-1"!), 1 < |a| < 4. 

They proved the existence and uniqueness of the constant mean curvature 
foliation {Mr} for SAF manifolds, where the mean curvature of Mr is (2/r) — 
{Am/r'^) + 0{r^'^). They showed that the geometric center defined as follows 
converges: 

C^y=lim^f^- , a = 1,2,3, (1.7) 

where z is the position vector of Mr in M. This definition is also motivated 
by the spatial Schwarzschild manifold in which the constant mean curvature 
foliation is {\x — p\ = r}. Another definition by Corvino and Schoen in 
[CS06| is defined for AF-RT manifolds by 



CrQ = lim / x'^iqui — qu j)vid(Ta 



{hiaK - hiii'g)dag 



\x\ ^ 



where we recall that hij = gij — 6ij , see also jB087[ IBLP03j . One applica- 
tion of this definition in jCS06j is the gluing theorem which allows Corvino 
and Schoen to approximate AF manifolds by solutions which agree with 
the original data inside a given region and are identical to a suitable Kerr 
solution outside a compact set. Their definition is a local-coordinate expres- 
sion which is natural from the Hamiltonian formulation and convenient for 
calculation purposes, but may obscure interesting geometry. 
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The main purpose of this paper is to prove that the intrinsic definition 
()1.4p is equivalent to the Corvino-Schoen definition (|1.8p . Moreover, for SAF 
manifolds in which the unique foliation of constant mean curvature surfaces 
exists, our intrinsic definition is equal to the Huisken-Yau definition (jl.7p . 
In other words, the intrinsic definition is a coordinate-free expression of the 
Corvino-Schoen definition and it generalizes the Huisken-Yau definition. 

Theorem 1. Assume {M,g,K) is AF-RT H^) . Then 

Ci = Ccs- 

Theorem 2. Assume (M, g, K) is SAF < f77g|) . Then 

Ci = Chy- 

We would like to make a note that Corvino and Wu |CW08j recently have 
a result about the equivalence of these definitions under the assumption that 
the metric g is conformally flat at inflnity with vanishing scalar curvature. In 
that special case, they are able to derive some explicit estimates. However, 
it seems that their approach cannot be generalized to AF-RT metrics. 

This article is organized as follows. In section 3, we prove a density 
theorem (theorem 12. 2p for AF-RT manifolds. The theorem is crucial for 
most of the arguments in the paper and may be of independent interest. In 
section 4, we discuss natural properties of the intrinsic deflnition. Theorem 
[Hand Theorem [2] are proved in sections 5 and 6, respectively. 

2. The Density Theorem 

Let (M, g, K) be a vacuum initial data set. We introduce the momentum 
tensor 

^ij ^ K'^ -Trg{K)g'^ . 
The constraint equations (II. ip then take the form 

divg(^) = 0, (2.1) 

and we define 

$(<7, vr) = + \{TrgTTf - \tt\1 , div,(7r)) . 

In the case that {M,g,TT) is AF-RT, Corvino and Schoen |CS06j prove that 
AF manifolds can be approximated by AF-RT manifolds in some weighted 
Sobolev space. Moreover, instead of requiring smooth solutions to the con- 
straint equations ()2.ip . their theorem works for solutions with weak regu- 
larity. Before we state the theorem, we need the following definitions. 
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Definition 2.1 (Linear and Angular Momentum). The linear momentum 
{Pi, P2, P3) and the angular momentum (Ji, J2, J3) are defined as follows. 

1 f 

Pj = — lim / TTi'jV'r, d(T„ 



1 



= lim / TTjkZ^ .iy^ dag, 

where give the rotation fields in M^, for example = ^^^^l^- — x^-^- 

Remark. The linear momentum is well-defined for AF manifolds and the 
angular momentum is well-defined for AF-RT manifolds. 

Definition 2.2 (Weighted Sobolev Spaces). For a non-negative integer k, 
a non-negative real number p, and a real number 5, we say f £ w'^f{M) if 



WlfiM) - 



E (\DV\P^''^^'Y daA <oo 



\a\<k 



where a is a multi-index and p is a continuous function with p = \x\ on 

M\Br,. 

When p = CO, 

ll/llw''=-(M) = E e5ssup|Z)"/|/,H+'5. 

^ ' rz, M 

\a\<-k 

Definition 2.3 (Harmonic Asymptotics). (M, (7, vr) is said to have harmonic 
asymptotics if (M, g, vr) is AF and 

g = u^6, TT = u^{C5X) (2.2) 

outside a compact set for some u,X tending to 1,0 respectively, where for 
any metric g, CgX is the operator associated to the Lie derivative Lxg 
defined by 

CgX = Lxg - divg{X)g. 

By the constraint equations ()2.ip , u and X in definition 12.31 satisfy the 
following equations outside the compact set, 

SA^^x = ( - |£5X|2 + ]^{Trs{CsX)fy, 

AsX' + Au-^jiCsX)] - 2u-^u'Trs{CsX) = 0, 

where Ui = = Asymptotic flatness requires u,X tend to 1,0, re- 

spectively, at some decay rate. Using the decay conditions on u and X, we 
have Asu = 0{\x\-'^-'^^) and AgXi = 0{\x\-'^-'^^). As shown in [B86], the 
asymptotic behavior implies that outside a compact set, 

u = l + -^ + Oi\x\-'~^), X' = ^ + Oi\x\-'-^), (2.3) 
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for some constants a and 6*. Note that this clearly implies that (M, g,7r) 
is AF-RT. Furthermore, if (M, g, it) is AF-RT, u"'^'^ and {X'y^ satisfy equa- 
tions with better decay: Asu°'^'^ = 0{\x\-^~^^) and AsiX')"'^'^ = Odxl-^-^''). 
Hence 

u'^'i'^ix) = ^ + 0(|x|-2-^), iXr'"'ix) = ^ + 0(1x1-2-^), 

(u'^dd^^^_^^ =0{\x\~^-'), (^(xr''''{x)-^^^ =0(|x|-3-^), 

(2.4) 

for some vectors c, which are quantities corresponding to the center of 
mass and angular momentum of {g, n). Since we assume that g and K satisfy 
the pointwise regularity, the above identities hold pointwisely. However, we 
can generalize the above discussions to the setting of the weighted Sobolev 
spaces and have more general results as follows. 

Theorem 2.1. [CSOHl Theorem 1] Let {gij-6ij,7rij) e W^f{M)xW^f_^{M) 
be a vacuum initial data set, where 6 G 1) and p > |. Given any e > 0, 
there exist ko > and a sequence of solutions {gk,T^k) with harmonic asymp- 
totics satisfying \2.2\) . 112. 3|] so that 

Moreover, the mass and the linear momentum of {gk,T^k) O'l"^ within e of 
those of {g, it). 

The theorem says that the solutions with harmonic asymptotics (|2.2p 
are dense among general solutions. More remarkably, the mass and the 
linear momentum which can be explicitly expressed for solutions with har- 
monic asymptotics converge to the original initial data set in these weighted 
Sobolev spaces. However, in the above theorem the center of mass does not 
seem to converge, neither is the center of mass Cj defined generally for AF 
manifolds. Therefore, we would like to modify their theorem and prove, in 
some weighted Sobolev space, solutions with harmonic asymptotics form a 
dense subset inside AF-RT solutions so that the centers of mass and the 
angular momentum converge. The precise statement is as follows: 

Theorem 2.2 (Density Theorem). Let{g-6,Tr) G W^f{M)xW^f_g{M) be 
a vacuum initial data set and (godd ^ ^even^^ ^ W^f_g{M\BR^) x w]_'^^s{M\ 
Bjig), where 6 G (^, 1) and p > 3. Given any e > and Sq G (0,(5), there 
exist R, ko = kQ^R), and a sequence of solutions {gk,T^k) with harmonic 
asymptotics satisfying i2.2\) . i2.3\) . ^2.4\ ) so that {gk,T^k) is within an e- 
neighborhood of (g, vr) in the WZf{M) X Wj^_g{M) norm and 

\\9k^'^\\w'^f_^ {M\Br) - ^' lkr''"llH/^'|_^- (M\Bfl) - ^' ^^'^ ^ - ^0- 
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Moreover, the mass, the linear momentum, the center of mass, and the 
angular momentum of {gk-,T^k) o,re within e of those of ((7, tt). 

We first briefly describe Corvino and Schoen's construction of the approx- 
imating solutions {gk,T^k)- 



Sketch of the proof of Theorem [2A[ Let {gk-,T^k) be 2-tensors cut off from 
the original solutions {g, n), gk = ikQ + (1 - ik)5 = S + ^kh, frk = CfeTr where 
^k is a smooth cut-off function 

{1 when |x| < k 

between and 1 when A; < |x| < 2A; 
when > 2k. 

is chosen so that \D^k\ 1^ c/k and ^ c/A;^ for some constant c 

independent of k. Then we let 

9k = ulgk, 

^k = Uki^^k + ^gk^k)- 

To simplify notation, we denote Cg^. by £, and we also drop the subindex 
k when it is clear from context. In order to find u near 1 and X near 
at infinity so that {g, vf) is a vacuum initial data set, we need to solve the 
following systems for u and X from the constraint equations (j2.ip 

fl = u-^(^- 8AgU + (^R{g) - Ivr + + ^ {Trgin + CX))^^^ = 

(diVg(7f))^ = U-^(^{diYg{lT + CX))^ + Au-^g^''uj{TT + CX)ik 

-2u~^UiTr^{T: + CX)^ i = 1,2,2,. (2.5) 

Consider the map T : (l,0)+VF^f(M)xVF!f (M) ^ VFlf_5(M)xVF°'2%(M) 

defined by T{u,X) = (/i, divg(7f)) . It is known that Z)T(i o) is a Fredholm 
operator of index 0. For the Fredholm operator of index 0, the operator is 
injective if and only it is surjective. However, it is not clear whether LT^j^ g) 
is surjective. Corvino and Schoen enlarge the domain and utilize initial data 
sets of the form {u^g + h, u'^{tt + CX) + with 

$ {u^g + h, u^{t^ + CX) +q) = (0, 0), 

where h and q are symmetric (0, 2)-tensors with compact supports. Then 
they prove that the operator L'<I>(g^^) maps surjectively onto VF_'f_^(-M) x 
VF°f_5(M) for p > 1 and 5 G (0, 1). 
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Since -DT^^ q) is Fredholm of index 0, we have 

W^fiM) X W^f{M) = Ker (DT(i^o)) ® Wi, 

W'j_s{M) X W^°f„5(M) = Range {DT^i,o)) ® span{yi, . . . , IV} 

where Wi is an A^-dimensional hnear subspace and {Vi, • • • , Vn} is a basis 
for the cokernel of DT^^i ^y Because L><&(-^ .^) is surjective, we can choose 
{ (/ii, . . . , (/iTV, (ZAf)} so that D^(^g^^-^{hi,qi) = Vi. Supports of those 
(hijQi) may not be compact, but there exist {hi,qi) with compact supports 
close enough to {hi,qi) so that Vi = D^(^g .f^-^{hi,qi) still span a complement- 
ing subspace for Range (DT(i q))- 

Let W2 = span{(/ii, gi), . . . , {h^, qw)}- W = Wi x W2 is a Banach space 
inside W^f{M) x W'^f{M) x W^f{M) x W^lf.^lM). Define the map T 
from ((1,0), (0,0)) +W1XW2 to VF1'2%(M) x VF°'2%(M) by 

T((u, X), (/i, g)) = <^{u^g + /i, u2(^ + CX) + g). 

-DT^^i o),(o,o)) is an isomorphism by construction. Hence the inverse func- 
tion theorem asserts that T is an isomorphism from a fixed (independent 
of k) neighborhood of ((1,0), (0,0)) to a fixed neighborhood of $(5, vr). Be- 
cause (0, 0) is contained in the image when k large, there exists a unique 
{{u, X), {h, q)) within that fixed neighborhood of ((1, 0), (0, 0)) such that 
$ {u'^g + h, n2(7r + CX) + q) = (0, 0) for k large. □ 

Note that the supports of h, q in the proof may not be uniformly bounded 
in k, but it is important in the proof of Theorem 12.21 that h, q have compact 
supports uniformly bounded in k. Therefore, we need to carefully choose 
the cokernel of DT^^i^q^ for k large. This choice is described by the following 
lemma. 

Lemma 2.4. V and W are Banach spaces. Assume that : V ^ W is 
a sequence of Fredholm operators and converges (in the operator norm) 
to some Fredholm operator S' . If W = RangeS" © W for some finitely 
dimensional closed subspace W , we can choose the cokernel of inside the 
cokernel W of S' for some C W , for k large. 

Proof. Since S' is Fredholm, there exists V' such that V = KerS" © V' . 
Consider for any bounded linear operator S -.V ^ W, the map 

Ts -.V ew' 

given by ts{v, w) = Sv + w. Then ts,, is an isomorphism for k large since Tg' 
is an isomorphism by construction and isomorphism is an open condition in 
the space of linear operators. We then have 

W = Sk{V') © W'. 
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Therefore, for any v £ V, Sk{v) can be decomposed uniquely into Sk{v') + 
Ski^) — Sk{v') for some v' £ V and Sk{v) — Sk{v') G W . Hence v can be 
decomposed uniquely into v = v' + {v — v') as well, where Sk{v — v') S W . 
Let Uk = {u £ V : Sk{u) S W'}. It is easy to see that Uk is a closed space 
in V and 

V = V'(B Uk. 

Note that KerSk is a finite dimensional subspace in C/^, so we can write 

V = V' ® KerSk Zk 

for some closed subspace C Uk- Sk{Zk) is closed in W, so there is 
Wk C W such that Sk{Zk) eWk = W' and hence 

W = SkiV) e SkiZk) (BWk = RangeSfc Wk. 

□ 

Corollary 2.5. The supports of h and q in the proof of Theorem \2.1\ can he 
chosen to he uniformly hounded for k large. 

Next, the key lemma (Lemma I2.8P used to prove Theorem 12.21 is an a 
priori type estimate for 2nd order elliptic equations Pv = f which have a 
symmetric property at infinity. Roughly speaking, if we know P and / are 
even (or odd, respectively) then we hope solutions v will be even (or odd, 
respectively) as well. However, this is not true generally because boundary 
values can affect solutions dramatically. For example, consider two harmonic 
functions ui , U2 in M.^ \ Bf{^ , 

1 C ' X 

ui = - — r and U2 = — rr- 

\x\ \x\-^ 

Both I nil and \u2\ tend to zero at infinity. However, ui is even and U2 is 
odd. They are solutions for different boundary values on the inner boundary. 
Nevertheless, in the case that the boundary value is very small, we will show 
the symmetry of the solutions is not affected much in the region away from 
the boundary. Before we state the lemma, we will give some definitions of 
the operators we consider. This class of operators is discussed in detail in 
[B86] . 

Definition 2.6. Let P defined as Pu = a^^{x)dfjU + U{x)diU + c{x)u he an 
elliptic operator. 

(1) P is said to he asymptotic (at rate t) to an elliptic operator P, 
Pu = a'^{x)dfjU + b'^{x)diU + c{x)u, if there exist q G (3, oo), r > 0, 
and a constant C such that over the region M \ B^, 

h'^ -a'^Ww'''^ + W -b'\\wO,g + \\c-c\\ o,,/2 < C. 

~T — 1 — T ^* —2 — T 
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(2) p°'^'^ is the odd part of the operator, defined on Ad \ Bji^ by 

Remark. We assume that P : W^^ — )• W^'^^ for s = 2 or 3, 1 < p < q and 
a non-integer a > 0, and define 

\\P\\op;M\Br = sup{||Puj||^^.-2,p(^) : llu^llw'-.P(M) = l,suppu' CM\Bb}. 

From Definition 12.61 we have 



\\Pw Pw\\y^,s-2,p^j^^ <\\a'^ "-'^^^wl-^iMXEn) 

+ 11^' - b'\\w°J?_^{M\Bn) + 11^ ~ ^K°J^^I{M\Br) ^ 
as i? — )• cxD. Therefore, \\P — P\\op;M\Bii — )■ as i? — )• 0. 

Definition 2.7. We say that a sequence of elliptic operators P^, PkU = 
a^^^^{x)dfjU + b^^^^{x)diU + C(^i^^{x)u, is asymptotic to P uniformly if given e, 
there exist R and ko such that 

ll"(fc) - "''^\\wb?{M\Bii) + ll^(fc) - ^'\\w1-^_JM\Br) 

Lemma 2.8. Let A5 he the Euclidean Laplacian in M^. Let {Pk} be a 
sequence of 2nd order elliptic operators asymptotic to A5 uniformly, and for 
s = 2 or 3, 1 < p < q, and a non-integer number a > 0, 

Pk : W'jliM) ^ VF^^'f (^)- 
Assume {vk} C Wff(M), {fk} C W^'^^{M) are sequences of functions. 
We also assume {vkY'^'^ G VF!f_jM \ S^jJ, \\{vkY'^'^\\wl:l{M\Bn,) ^ 

— )• 00, and {vky^ satisfy Pk{vk)°'^'^ = fk- Then there exist R large and ko 
large such that 

wherec = c{s,p,a, \\Pk-As\\op;M\BRj,R = R{s,p,a, \\Pk-I^5\\op;M\BR^) ,ko ■ 
ko{R). 



Proof. Because {vk)°'^ is not a global function, we multiply it by a smooth 
function 0, 

{0 when \x\ < R 

between and 1 when R < \x\ < 2R 
1 when |x| > 2R, 

where (p satisfies \D(p\ < c/R, \D'^(p\ < c/R'^ for some constant c independent 
of R, where the real number R will be chosen later. We have 4'{vk)°'^'^ G 



On the Center of Mass of Isolated Systems with General Asymptotics 11 

W"!f_a(M) and 

\odd\ _,f, .ij I y -.odd 



By |B861 Theorem 1.7], there exists a constant ci = ci(s,p, o) such that 

U{v,r'^\\w-f_^iM) < ci||A5(<^(r;fc)°'^'^)||^.-..,(^,). 

We then consider the difference between As and Pk by viewing them as the 
operators from H^l'f„^(M) to Wl3^'^(M). By the remark after Definition 
2.61 the operator norm — ^5\\op;M\Bii tends to zero as i? — )• cxd uniformly 
in k because Pk is asymptotic to uniformly. Therefore, 



< 



Ci\\Pk - As\\op-,M\BjHvkr^^\\wlf_JM)+Cl\\Pk{Hvkr^^^^^ 



Choose R such that ci||Pfc — As\\op-m\Br ^ ^ ^^'^ absorb that term to the 
left hand side. Then 

mvkr''''\\wif_jM) <2ci\\Pk{c^{vkr'''')\\wi-i':{M) 

2 \odd\ 



<C2\\fk\\w^-2,PjM\B,,) + C2\\{D'^){vkm\wlf_t{Aa) 



<C2\\fk\\wif_f(A4\BR) + <^3\\{Vk)°'^'^\\ws_-l,P^AR)^ 

where Ar = {x : R < \x\ < 2R}. Because \\{vk)°'^'^\\^rs-i,P(^j^'^^^ ^ — )• as 
/c — >• DO, for e = 1/{2R), there exists ko such that for all k > ko 

ll(^^fc)°''llvF!7!-(A,) < 2ii||(^;,)'"^'^||^.-i,p(^^) < 2Re = 1. 

As a result, we have the Schauder-type estimate 

\\M°'^'^\\wt'^^JM\2R) < C2||/fc|lvi/»-2.P(jvf\B^) +C3 for all k > ko. 

□ 



Proof of Theorem \2.Si 

1. Estimates on u°'^'^ and [X^y^"^. We construct {gk,T^k) as in Theorem 12.11 
in the form 

9k = Uk9k + hk, vffc = ul^TTk + Cg^Xk) + qk- 

Recall that k is the radius of which we cut off the original data. Again, we 
drop the subindex k when it is clear from context. By Theorem l2.lt u and 
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X exist and satisfy the system of the constraint equations ()2.5p . From the 
constraint equation (12. Sp for u in M \ Br^ , we have 



= AgU-l(R{g)-\7r + CX\l + ^{T^7r + CxY)u 



dx^dx^ ' ^ ^ ("c?x*^ ^ dx^ 
-i (ii(ff) - Ivr + + i (Tr^(7r + " 

= Piu. 

On M \ , n°'^'^ satisfies the following equation, 
P^u°'^'^ = {p^u){x)-{Piu){-x)-Pf'^{u{-x)) 



odd 



R{g)-\7T + CX\l + ^{Trg{7r + CX)] u{-x) 

{gHx)-9H-x)) Q^<-^) 
ll ^{x)[^^g'^{x)^/l{x)^ 

-1. 5 s /T. ■.\\ d 

1 1 ( 



/1+/2, 



where 



1 / -1 \ odd 



h = ^[-\T^ + CX\l + -TTg{TT + CX)j U{-X), f2 = PlU°''''-h. 

fi contains the terms involving X. We will use a bootstrap argument to 
improve its decay rate. /2 contains the terms which have expected good 
decay rate already, such direct calculation tells us 

We emphasize that through out this proof, c is a constant independent of k. 
From the constraint equations (|2.5p for X, we have 



{diVgiCX))^ + diYg{7T), + 4u-'g^^uj{iT + CX). 



ik 



- 2u~^niTr§(7i- + CX) = 0. 
If we compute the first term in local coordinates, we have 
(div^(£X)) . = (AgX)^ + RicikX'' 
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We define P2Xi = g^^ gSdx^ ^^"^ then P2Xi = Fi, where Fi contains the 
remainder terms from above identities, 

-div^(7i-), - iu-'^g''^Uj{TT - CX).^ + 2u~^UiTx-g{T^ + CX). 
Then we have 



= ai + hi, 

where 



ai = - {di\g{T^)i)°'^ , bi = P2 {Xi)°'^'^ - ai, 

where bi contains u°'^'^ and X"'^'^ and we will bootstrap to improve its decay 
rate. A straightforward calculation gives us 

M W°_f_,{M\BR^ ) ^ 1 1 6i 1 1 ^o,p_^^ (j^^^^^^ ) < C. 

From the above, we derive the system 

Pin''"" = /i + /2, 

P^iXiy^"^ = ai + bi, i = 1,2,3. 

We apply Lemma 12.81 (with a = 26 — 1 < 6) for each equation because Pi 
and P2 are obviously asymptotic to and 



1 1 (-^i ) fc"''^ 1 1 w^-p (m\Bh„ ) - I K-'^i ) fc'^'^ 1 1 {M\BRg ) < 2 1 1 (Xi ) 1 1 ^2,p (^j) ^ as k 
Because — 1 — a = 25, by Lemma 12.81 there exist Ri and ki such that for all 



k > ki. 



1 1 (^fc 1 1 w^fg (Af J - 1 1 ) "k^"^ 1 1 W^fg (A/ ) - 



Once we derive these estimates, the decay rates for fi and bi are improved. 
The bootstrap argument allows us to conclude that for some R2 > Ri, 
k2> ki, 

||(nfc)°'^'^||^3,._^(^,^^^^) < C, \\iX^)n\^^_._^^M\B,^) < C for all k > k2. 

Therefore, for any given e and 6q G (0,(5), there exist R and /cq so that for 
all k > ko. 

Furthermore, by Corollarv l2.51 the supports of (/i^, q^) are uniformly bounded 
in k. Hence we have 

1 1 7,odd II ^ 

Similarly, we have the estimate for f . 
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2. Convergence of the Center of Mass and the Angular Momentum. To 

prove the center of mass and the angular momentum of (5, vf ) converge to 
those of {g,7r), the same idea of proving convergence of the mass and the 
linear momentum in [CS06j is employed. 



\Cf{g)-Cn-9)\< 



+ 



+ 



Cf{g) 



1 



(R 



Rg9ij)y{a)^l dag 



cm 



1 



, R 



-Rgij)Yl^^vlda-g 



The first and the third terms can be written as integrals over > r} by 
the divergence theorem. We will also use the fact that > r} is centrically 
symmetric, i.e. x, —x € {|x| > r}, to estimate those integrals over {|x| > r}. 
For the first integral, 



Cf{g) 



1 



\x\=r 



Ri3 - i;R99ij)Ylc,)^idag 



x\>r 



X 



"Rg-2x''{g'^ -6'^){Rj--Rggij) 



+ {R^- Rgg,j)Yi^^Tig^^ 



dvol 



a- 



(2.6) 



Let H = {x £ M \ Br : x^ > 0} be the half space. Then we can rewrite the 
above integral as follows: 



Jh 



H 



2x^{g^^r''\R,--Rgg,,)^dx 



L 



noc 

R,j 



odd 1 



^Ro/^g,^--Rg{-x)gff^)^ 



1 



+ - 2^g9ij)V9' 



■odd 



dx 



+ 



YL^ R 



H 



■yodd 



\Rf'g^J-\RgC^^f)v%^dx 



+ Yi^^{R,{-x) - -Rgi-x)g,,{-x)) {Tl^r^^^ + rl^{-x)^' 



■odd 



H 



dx. 



We substitute Rg in the first integral by using the constraint equation Rg = 
— |(Trg7r)^ + \TT\g. We then bound above integrals symbolically and use the 
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Holder inequality: 

ci / \x\(\T:'-'^M + \ii\''\g°''''\ + \g"''''\\D''g\ + \^ 

J H 

+ cij^ \x\\\D\g"'"')\\Dg\ + \D^ g\\D{g"'"')\ + \D^g\\DgU'"'\) dx 



^ n I Uneven 1 1 ll^ll i ll„l|2 \\ „odd\ 



/2,p 



—1—6 —6 / 



where ci and C2 are constants independent of g, vr, and r. The weighted 
Sobolev norms above are over the region {|x| > r}. Similarly, 

cm - [ {R.j-lR-gij)Yi^)ulda-, 

J\x\=r ^ 
^ „ M I - even ii ||-|| ill;s^l|2 || -odd 1 1 

T.odd\\ ii-ii \„l-2(5 



nr^'i^^. mw-Ar' . (2-8) 

— 1 — 6 —5 / 

For the surface integral, we can assume k ^ r (recall k is the radius of 
which we cut off the original data) and then g = u^g on {|x| < 2r}. Hence 
we have 



/ [Rii - lRa9tj)Y(a)^^g dag - / {R^j - ]-Rgij)Yl^^vl da-g 

J\x\=r ^ J\x\=r ^ 

= 1^1^^ ((1 - y'')R^J - ^(1 - u')R99^J)Yi^)'y'g dag 

+ ^ ^ (u^ (Rij - Rij )-^u''{Rg- R)gij ) Y^^^ dag 
< C6max(|l -nllx]-^) max (\D^g\\x\^+^)r^'^^ 

\x\=r V / |x|=r V / 

+ C6 max _ |^|2+5^2-5_ 

\x\=r 

Recall the Sobolev inequality for weighted Sobolev spaces (see [B86]). For 
n — sp < where n = dimM = 3 , we have ||?i;||^^o,ooj.^^^^ < c||w|||ys,p(7;//). 

Therefore, for p > 3, 

mayi\D^{g -g)\\x\'^+^ < max \D^{g - g)\\x\'^+^ 

\x\=r M\Br 

= \\D^{9 - 9)\\w'j^_,{M\B^) ^ cp'(5 - 9)\\w^_:^,_,{M\Br) 
< 49 - 9\\w^_,V(^M\B,y 
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Similarly, using the Sobolev inequality on the other term, we derive 



/ {Rij - lRg9tj)Y(a)^^g dag - / {Rij - ]-Rgij)Yl^^ul da-g 

J\x\=r ^ J\x\=T ^ 

< C7(||l - ^^IIh/!-^(m\B,) + \\9 - 9\\w^^,P(^M\Br))''^'^- (2.9) 

To conclude the argument, for given e, let Ri be a constant so that the right 
hand side of (j2.7p is bounded by e/16. For this e, there exist R2 and k > koso 
that \\1-Uk\\^^3,pf^j^jy \\g-gk\\^^3,p^j^^y IK-7rfc||^i,p^^^^^p \\'^k"'''\\w^-^_g{M\BR^) 
and ||5fc'^'^||yj/3,p (M\Bji ) small, for all k > ko. Using estimates (12. 7[) . 
()2.8p . and (|2.9p . we can conclude that, for all k > ko, 

\Cf{g)-Cf{-g,)\<e. 
To prove that angular momentum of g is close to that of gf, we notice that 



J„-/ TTjkZl^^iy^ dag = diVg{TTjkZl^^) dVOlg 

J \x\=r J \x\>r 



|x|>r 

The first term is zero since TTjk is a symmetric tensor and Z'^ is a Killing 
vector field. The other terms after integration are bounded by Cr^^^^. Then 
the rest of the argument works the same as the case of center of mass. □ 



In the proof of convergence of the center of mass and the angular momen- 
tum, we showed that the limits at infinity can be approximated by surface 
integrals at the finite radius and the difference for surface integrals at the fi- 
nite radius is arbitrarily small. Since we will use this argument several times 
through this paper, we formulate the argument into the following lemma. 

Lemma 2.9. Let F{g) he a vector field depending on x , gij , D gij , D"^ gij 
smoothly defined in M \ i?_R(, . Let a,b > 0, and let r be the radius. As- 
sume that the following estimates hold for k > r, and for g,gk satisfying the 
assumptions in Theorem \2.2\ (Density Theorem). 
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(1) 



x\>r 



divgF{g) dvolg 



< c IItt' 



+ P" 



I n"'^'^ I 

W^f_,(M\Br)^\9 \\w'^f_^{M\Br) 



+ \\9°'^\\w^_'P_^{M\Br)\\S - 



W^f{M\Br) + llf ^\^W^'^{M\Br) 



(2.10) 



(2) 



J\x\=r J\x\=T 



9k '^'^9k 



< c| 111 - u\\w'^jV^^M\Br) + \\9 - 5'fellH/^'f_^(M\B,) 



(2.11) 



Then given e > 0, there exists ko such that 

/ F{g) ■ Vg dag = lim / F{gk) ■ Vg^^ dag^ + e, for all k > ko- 

J\x\=r r^°°J\x\=r 



lim 



Remark. If we replace the normal vectors and the volume forms of the 
integrals in the above assumptions by those with respect to the induced metric 
in the Euclidean space, the analogous result is the following: assume the 
following estimates hold 



(1) 



|3;|>r 



divsF^g) dvolo 



+ M 



odd I 



W^_-[_,{M\Bry\3 llw!f_,{M\B.) 



+ \\9°'^\\w^'P_^{M\Br)\\S ^W^f{M\Br 



+ 115-^11^2,. 



W^'P{M\Br) 



(2.12) 



(2) 



|j;|=r 



F{g) ■ vq dao - / F{gk) ■ vq dao 



(2.13) 



Then given e > 0, there exists ko such that 
lim / F(g) ■ vq dao = li™ / P{9k) ' dao + for all k > kg. 



18 



LAN-HSUAN HUANG 



3. Properties of the Intrinsics Definition 



Bartnik |B86l Section 4] proved that the ADM mass is a geometric in- 
variant and satisfies natural properties. Along the same lines , we will show 
that the intrinsic definition of center of mass ()1.4|) is well-defined and has 
corresponding change of coordinate under the transformation at infinity. We 
will first show that the intrinsic definition is robust in the sense that we can 
integrate over a general class of surfaces, and those integrals converge to the 
same vector. 

Proposition 3.1. Suppose {M,g,K) is AF-RT fTgj. Let {Z^fcl^i C M 
he closed sets such that the sets Sk = dDk are connected two-dimensional 
-surf aces without boundary which satisfy 

rk = inf{|x| : x £ Sk} — )• oo as k ^ oo, (3.1) 

r^^area(S'fc) is bounded as k ^ oo, (3.2) 

vol{Dk\D-} = 0{rl), 

where = Dk H {—Dk} and 3~ is a number less than 3 . 

(3.3) 

Then the center of mass defined by 
is independent of the sequence {S^}- 

Remark. The first two conditions (j3.ip , (j3.2p on S^ are the conditions consid- 
ered by Bartnik [B86] to ensure the ADM mass is well-defined. The volume 
growth condition (j3.3p allows us to consider a general class of surfaces which 
are roughly symmetric; that is, the non-symmetric region Dk\D~^ of has 
the volume growth slightly less than the volume growth of arbitrary regions 
in M. 



Proof. As in (12. 6p . the divergence theorem gives us 

(i^^i - \Rg9ij)Yl^)vl dag = {R,j - ]^Rg9ij)Yl^)vl dag 

+ [ ^""Rg + 2x°(5^^- - 5^^){R,i - \Rgg^.j) + {R^, + \Rgg.,)Y^^^n^ dvolg. 

We can decompose the integral over Df^ \ Di = {Df^ \ D^} U {D^ \ D^} U 
{D^ \ Di} into three integrals, where = Di U {~Di}. The first integral 
over {Dk\D^} converges because of the volume growth condition (|3.3p . The 
second integral over {D^ \ Df} converges because the region is centrically 
symmetric and the initial data set is AF-RT. After taking limits, the right 
hand side has a limit independent of the the sequence {Sk}- 
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□ 



Assume that {x} and {y} are two AF coordinates on {M\Bjig , g). Assume 
that F is the transition function between these two coordinates and y = 
F{x). It is shown that the only possible coordinate changes for AF manifolds 
at infinity are rotation and translation |B86l Corollary 3.2]. More precisely, 
there is a rigid motion of R^, (Oj, a) € 0(3, R) x R so that 

\F{x) -{Ox + a)\e W^'^^iR' \Bii,). 

Similarly, if we use the same argument in that corollary for AF-RT mani- 
folds, we derive 

Because the center of mass is a quantity depending on the coordinates, 
we now show that centers of mass in {x} and {y} coordinates have the 
corresponding translation and rotation. An interesting phenomenon is that 
compared to rotation, translation is a more subtle rigid motion. If the 
translation a is not zero. Density Theorem (Theorem 1 2. 2 p is involved in the 
proof. 

Theorem 3.1. Let {x} and {y} be two distinct AF-RT \1.3\) coordinates 
for (M, g, K) satisfying the change of coordinates as we described above. 
Assume that Cj^x and Cj^y are the centers of mass defined by the intrinsic 
definition in these two coordinates, then 

Ci,y = OCj^x + a. 

Proof. The metric g in the cotangent spaces induced from {x} and {y} can 
be written locally as 

ds"^ = gki{x)dx^ dx'' = gij{y)dy'dy^ 

= gij{F{x))d{Oix'' + a')diOjx^ + a^) + 0{\x\"^). 

Expanding the above terms, we get 

gkiix) = gij{F{x))Oloj + ei{x) 

where the error term ei G VF^'^^ and e°'^'^ G W]_'^. Then a straightforward 
calculation gives us 

d'gki _ d^g,, 



dyPdyi 

and the formulas for Ricci and scalar curvatures, 

Rki{x) = Rij{F{x))0i0i + e3{x) 
R{x) = R{F{x)) +e4(x) 
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where eq{x) G 0{\x\-^) and ef'^{x) = 0{\x\-^) for q = 2,3,4. To calculate 
Cj^y, Propositfon 13.11 allows us to integrate over {\y — a\= r}, 

We would like to replace y = F{x) in the above identity. First we have 

|y|25"i - 2y"y^ = - 20'^x^Olx'')6'^' + 20x ■ a6'^' 

-2{a''Oix^ + O^x^a') + 0(1), 
- _ Ojx^ 



\y - a| 



+ 65, 



where 65 G and ef"^ G W_^. Using the facts that Rij{F),R{F) and 

g{F) are asymptotically even 

+ / (R^J{F{x)) - \R{F{x))g,,{F{x))) ( - 2a''Olx'^) dao 

J\x\=r ^ ^ ' ^ \^\ ' 

+ / {%,{F{x))-\R{F{x)Yg,,{F{x))\i20x-ab^'-20^^^^ 

J\x\=r ^ ^ ' ^ ' \^\ 

+ 0(r-i). 

We claim that the first integral 1\ converges to IGTrmO^Cj ^, the second 
integral I2 converges to IGvrma", and the third integral 13 converges to 0. 

h = \ - iii(F(x))g,,(F(x)))0^0/(|x|20^ -20^xV)f^dao 

J\x\=r ^ ^ ' ^ ^1^1 

= 0^1 _ (i?fc,(x)-^/?,(:E)<7fc/(^))(kl'5'''-2xV)^(iao+0(r-i) 

= 167rmO^Cf^^ + 0(r-i). 
By using the formula for ADM mass (jl.Sp . 

r 1 ' 

I2 = a" ^ ^ - -Rg{x)gu{x)) ( - 2x'=^) do^ + O(r-i) 

= levrma" + 0(r~^). 
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To see that I3 converges to 0, we simplify the expression by letting b = O^a, 
and then 



1^ = 0^- 



Rm{x) - ^Rg{x)gM{x)) (2x ■ bS^^ - 2x''5'=) ^ da^. 



By the fact that (2x ■ b6^^ — 2x^6^) ■ ^ = 0, the scalar curvature term 
vanishes after taking a limit, so we only need to prove 



[ Rh{x)(2x- 

J \x\=r ^ 

It is sufficient to prove 

Rki{x)(2x 



2x^b^]^dao = 0{r-^). 
' \x\ 



2xH^y dog = O(r-i). 



A straightforward calculation shows this identity holds if the metric is con- 
formally flat outside a compact set, so we would like to approximate g by 
g which have harmonic asymptotics, and then apply Lemma l2.9i We need 
to check if conditions ()2.10p and (j2.1ip hold. Using the divergence theorem, 
we have 



divn 



|x|>r 



Rki{x){2x-b6'^'' -2x'^b'') 



dvolr. 



|a;|>r 



{dwgRki){2x ■ b5^^ - 2x^b^) + Rki{2b^5^^ - 25^^b^) dvolg 



sum to =0 



Rki[2x 



2xmg'^TZnd^o\g 



|a;|>r 



Using the second Bianchi identity to the first term on the right hand side, 
symbolically, the above integral is bounded by 



div„ 



|a;|>r 



Rki{x){2x-b5P^ -2x^b^) 



dvoL 



Hence condition (|2.10p holds. Condition (j2.1ip holds because the Sobolev 
inequality implies 



/ Rki{x) (2x ■ b5^'' - 2x^b^'] vi dog - [ Rki{x) (2x ■ b5^'' - 2x' 

J\x\=r ^ ' J\x\=r ^ 



\wt:l(.M\Br) 



\9 - 9\\w'^>v 



w^!:l{M\Br) 
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By Lemma 12.9^ we conclude 

lim / Rkiix)(2x-b5^'' -2x^b'')iyida„ 
= lim / Rki (x) f 2x • b6^'' - 2x^b'') ui da-„ + e = e. 

Since e can be chosen arbitrarily small, converges to 0. Therefore, 



^ly = 1™ + ^2 + h) = O'^C^^ + a". 



□ 



4. The Corvino-Schoen Center of Mass 



Corvino and Schoen [CS06j defined the center of mass (jl.Sp for AF-RT 
manifolds. In this section, we will show that the intrinsic definition we 
propose is actually equal to the Corvino-Schoen definition (Theorem [T]). In 
other words, the intrinsic definition (|1.4p is a coordinate- free expression of 
the Corvino-Schoen definition. 

Assume g is the approximating solution in Theorem 12.21 (Density Theo- 
rem). A straightforward calculation shows Ci{g) = Ccsid)- As proven in 
Theorem [2]2l Ci{g) can be approximated by Ci{g). If we prove that Ccsio) 
can also be approximated by Ccsio), then Theorem [1] follows. 



Proof of Theorem 1. First notice that the Corvino-Schoen definition is equal 
to the following expression where the normal vector and the volume form 
are with respect to the induced metric in the Euclidean space 



1 



^lim^ j x" Y^igij,i - gii,j)vi da[ 

J\x\=r 
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A direct calculation shows 
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■J\x\=r2 J\x\=r2 ■ 

/ a;" - 9ii,j)'^l dao- ^(^^ 

•^|a;|=r-i -J J\x\=ri • 



ri<\x\<r2 



- 9ii,jj) dvolo 



It is easy to see that the conditions (j2.12p and (j2.13p in Lemma 12.91 (and 
the remark afterward) hold because 'Ylij{9ij,i3 ~ 9ii,jj) the leading order 
term of Rg, and by the constraint equation (j2.ip . 



/ ^"Yl (9ij,ij - 9ii,jj ) c^volo 

J\x\>r ^ 



<^ „ / ||„|| ||„eiien 



+ 



\9°'^'^\\w^j^_^(M\Br)\\9 ^\\w^j^{M\Br)) ^ 
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and 



( / x°'^{9ij,i- guj^ldao- j '^{hiai'o 

\J\x\=r J\x\=r j 



hiiUo) dao 



\x\=r 



H],i- 9ii,j)^odao- / ^{hiaT^o - hiiu^) dao 



\x\=r 



2-S 



< c\\9 9\\w':^f{M\Brf 

Then we can apply Lemma 12.91 and have 
1 



lim 
167rm r^co 

1 



J \ x\=T ■ ■ 



- 9ii,3)K d(^o 



lim 

IGvrm r-5>oo 

CSsi9) + e. 



/ ^" ^i9ij,i - 9ii,j)^o dcro 

J\x\=r -J 



'^{hiaJ^o - hiii'o)dao 

i 



+ e 



□ 



5. The Huisken-Yau Center of Mass 



In the case that {M,g,K) is SAF ([L6l), Huisken-Yau |HY96j and Ye 
|Y96] proved the existence and uniqueness of the constant mean curvature 
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foliation in the exterior region of M, if m > 0. Huisken and Yau used 
the volume preserving mean curvature flow to evolve each Euclidean sphere 
centered at the origin with a large radius and they showed that the Euclidean 
sphere converges to a surface with constant mean curvature. Furthermore, 
they proved those surfaces with constant mean curvature are approximately 
round, and their approximate centers converge as per equation (I1.7P : 



where {Mr} are leaves of the foliation and z is the position vector. They 
define this to be the center of mass for M. Instead of using volume preserving 
mean curvature flow, Ye used the implicit function theorem to find a surface 
with constant mean curvature which is perturbed in the normal direction 
away from the Euclidean sphere Sr{p) for some suitable center p and for the 
radius R large. Although it is not emphasized in his paper, p can represent 
the center for each constant mean curvature surface because this surface is 
roughly round. Furthermore, we will show that p converges to Ccs and 
use this fact to prove Chy = Ccs- Then Cj = Chy follows because the 
Corvino-Schoen center of mass Ccs is equal to the intrinsic definition Ci 
by Theorem [H Before we prove Theorem [21 we will need a technical lemma 
(Lemma 15.11) which suggests that the center of the surface p converges to 
Ccs- 

Let v be the outward unit normal vector field on Snip) with respect to 
the metric g. As R varies, v is well-defined in a tabular neighborhood 
of Sr{p). Therefore, the mean curvature at y G Sji{p) is div^^^^^z^ = 
diYgV, the divergence operator of the ambient manifold M. Recall that 



HY 



= lim 



a = 1,2,3 



(1.7) 




V\y -p\ 



V = 




yk _ pk Q 

\y - P\ dy^ 



+ 0{R-'^). 
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A straightforward calculation gives us the mean curvature on Sr{p) is equal 
to 

2 4m 6m(y — p) ■ p 9m? 
diVoZ/ = r - -: TTT H ■: n h 



.y^-p' Pii{y) 1 , .y^-p' , „ 

-PijAy)i r ~ I T + nPiid(y)-\ r + -^O' ^-^ 

\y-p\ \y-p\ 2 ly-yl 

where |£'o| < -^(l + \p\) for some constant c depending only on the metric 
9- 

Lemma 5.1. For R large, 



\y-p\=R 



[y -p )[ ^Pij,k \y-p\3 I ° 



J\y-p\=R V / 

+ / (y - P ) - Pij,ii r - 1 r + 7;Pii,ji r d(^o 

J\y-p\=R V \y-p\ 2 |y-p|y 

-8m7rCg5 + 0(ii"^), q = 1,2,3. (5.2) 



Proof of Lemma \5.1\ We denote the first integral by, for a = 1, 2, 3, 

In the proof, we will rewrite X°(i?), and then some cancellation allows us 
to rearrange the left hand side of (j5.2p so that it has an expression corre- 
sponding to C^g. 

Since the coordinate is only defined outside a compact set of M, we can use 
the divergence theorem only in the annular region A = {R < \y — p\ < Ri}, 



2 7aV \y-pr 



2 V \y-p\^ J,k 

~o / Pij^H 1 12 "^°^0 

2 Ja V 12/ -pr /,fc 

+ n / \ 12 ^^°l0- 

2 Ja V |y -pr y,i 
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Using the divergence theorem and simphfying the expression, we get an 
identity containing purely boundary terms 

X°(i?i) = X°(ii) + - J'^iR), for all Ri > R 

where 

2Jiy_pi=ji \y-p\ 

J\y-p\=R \y-pr 

1 /" y"' — p"' J 1 /" y^ — p^ J 
+7; / Pa-, r"^o + 7T / Pia-, r«o-o- 

2 J\y-p\=R \y-p\ 2 j\y-p\=R \y-p\ 

To prove that X°(i?) = J'"'{R), we would like to apply Lemma [2.91 It is easy 
to check that the conditions ()2.12p and ()2.13p hold, so we get, for a = 1, 2, 3, 

X"(i?) - J'^iR) = lim (X"(i?i) - = lim (X?(i?i) - J^{R^)) + e, 

Rl^oo 5-00 

where X? and denote the integrals that we obtain by replacing pij by pij 
in X and ^7°, where are 0(|y|^^)-terms in the approximating solutions gij 
in Theorem 12.21 (Density Theorem) . More precisely, gij has this expansion 



and Pij is defined by 



gij = + +0(|2/r'). 



Pij = gij - (1+ 



Because = |^ + 0(|y|~^) as in ([2^ . it is easy to see that 



lim {I^{R^)-Jp?{R,)) = 0. 

Therefore, we conclude X"(i?) = J"{R). We then replace X"(i?) by J'^iR) 
in the identity (|5.2p and derive that the left hand side is equal to 

iff y^ -p^ 

- o / (y" - p°')iPij,i - Pii,j)i r 

y^ — p^ y°' — p°' , \ 
Pia-, r ~ Piii r ""^o • 



y-pi 



=R \^ ~ P\ F ~ P\ 
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We rewrite the above integrals into the summation of the center of mass 
(jl.Sp and the remainder terms. Then exphcit calculations yield 



y3 — pj 

- P"){9ij,i - 9ii,j)i r dao 

\y-p\=R \y-P\ 



yl _ pi yO- _ pa 

9ia-. f — gu— — dao 



\y-p\=R \X — p\ 



\X — p\ 



+ 



iy" -p") {{9ij,i -Pij,i) - i9ii,j -Pii,j)) 1 r^^o 

\y-p\=R \y-p\ 



yi — pi ya — pa 

{gia - Pia)-. r - {gu - Pii)— — dao 



\y-p\=R \X-p\ 

8TTmCgs + 0{R-^). 



\x — p\ 



□ 



Proof of Theorem\^ Let -Fp,_R : <S'i(0) — t- M be an embedding defined by 
y = Fp,R{x) = Rx + p, that is, Fp^ji{Si{0)) = Sr{p), the Euclidean sphere 
centered at p with the radius R in M. We consider the perturbation along 
the normal direction on Sr{p) defined by S = {y + X(j)i> : y £ Sr{p)} for 
a parameter A > 0, and for (p £ C'^'"'iSji{p)) with ||(/)||c.2,q < 1. We denote 
the mean curvature on S by H{p,R,\(j)). Using this notation, H{p,R,0) = 
divgi^, the mean curvature of Sr{p). By Taylor's theorem for mappings 
between two Banach Spaces, 

H{p,R,X(j)) = H{p,R,0) +dH{p,R,0)X(p 

+ (1 - s) [(fH{p, R, s{X<p)) (A(/>, A(/>)) ds 

where dH and d^H are the first and second Frechet derivatives in the (p- 
component. In our case, dH{p, R, 0) is the linearized mean curvature oper- 
ator on Sji{p), i.e. 

dH{p, R, 0) = As^^p) + \A\l + Ric{i^, v), 

where Ag^(^p^ is the Laplacian operator on Sr{p) with respect to the induced 
metric from g, A is the second fundamental form on Sr{p) and Ric{-, •) is 
the Ricci curvature of M. Then 

H{p,R,X4>) = H{p,R,Q) + As^(^p)X4>+{\A\l + Ric{v,v))X(t> 

+ {l-s)(^d^H{p,R,s{X(l)))X'^(l)(f)^ds. (5.3) 
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In [HY96] ■ the estimates on the eigenvalues Ai and A2 of A are derived 
and 

\A\l = A? + Ai = + 0{R-'), Ric{v, v) = 0{R-^). 
For the second Frechet derivative in the Taylor expansion, we have 

d^H{p,R,s{X^))X^^^=^H{p,R,t{\ct>)) 

The right hand side is the second derivative of the mean curvature of the 
surface {y + s{Xcj))u : y S Sji{p)}. For R large, the unit outward normal 
vector field on {y + s(A(/>)z^ : y £ Sr{p)} is close to v, and a straightforward 
calculation gives us 



d' 



,2 



—H{p,R,t{X^)) 



<cX\\R,,ki\\AM' + \AM\D^cl,\ + \Af\cl,(^ 



where the constant c is independent of p, R, (j). Let G and Ei be defined as 
follows, where G is the lower order terms of the mean curvature of Sn^p) 
from (O), 



1 {y' - p'W - p>){y^ - p^) ^ ^ {y' - p')iy^ - p') 
Giy) = -Pij,k , To h2pij 13 

- p' Pa , 1 y^ -p' 
\y-p\ \y-p\ 2 \y-p\ 

and 

Ei{y) =Eo + {\A\l - (A0) + Ric{u, u){X(^) 

+ (1 - s) {(fH[p, R, sA0)) (A(/>, A(/>)) ds. 

G will give us Ccs as indicated in Lemma l5.ll and Ei{y) is an error term 
bounded as follows for some constant c independent of R, (p. 

C ,^ , C /,,,,, 9, , , ^2 \2| 



From identities (15.11) and (15.3 



H{p,R,X(l>) = -^-^+ \/ ^ + ^+G{y)+X^Sn{p)<P+J^^<P+Ei{y). 

To find the surface with constant mean curvature, we need to find ii, (p so 
that 

9 4?77 

H{p,R,XcP) = --^. 

It is equivalent to solving 

Q ^ My)-P + 9^ + c(^) ^ ^^^^^^^^ ^ + (5 4) 
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We pull back the equation (15. 4p via the map i^p,_R, and we get the following 
equation on 51(0), 

+Ei o Fp^nix) (5.5) 

for p,R,7p{x), where A^^^q) is the Laplacian on S'i(O) with respect to the 
pull back metric and tp = cp o Fp^R^x) = 4>{Rx + p) is the pull back of (/>. 
Define the operator L : C72'"(5i(0)) ^ (:70'"(5i(0)) by 

L = -Ao-2, 

where Aq is standard spherical Laplacian on Si (0) in . Because the metric 
g is asymptotically flat, the difference between Aq and A^^^q) is small and 
can be treated as the error term. Therefore, the identity (I5.5p is equal to 



= — ^ + + GoFp^R{x) - -^ALV'(x) + El o Fp^R{x), 

where Ei has the same bound as Ei. We let A = i?^" for some fixed 
a G (0, 1) and multiply on both sides of the above equation, then 

Furthermore, since D^ip = (Dy(j))R, 



In order to find p, R and ip to solve (j5.6|) . first we perturb p = p{R, ip) so 
that the right hand side of ()5.6p is inside RangeL for any R and tp. We will 
also show that p = Ccs + e where e is the error term containing lower order 
terms in R and Wi^W- Second, using an iteration process and the Schauder 
estimate, we can find a solution -p for R large. 

1. Perturb the center p. L has a kernel equal to spanjx^, x^, x^} be- 
cause translation preserves the mean curvature. Since L is self-adjoint, 
C°'"(5i(0)) has the L^-orthogonal decomposition C°'"(5i(0)) = RangeL 
span{x^, x^, x^}. We would like to find p so that the right hand side of (|5.6p 
is orthogonal to span{x^, x^, x^}. That is, we want to find p so that for 
a = 1,2,3, 

/6mx • p 9m^ \ 
+ [ x" (r^+'^G o Fp^Rix) + R^+'^Ei o Fp^R{x)) dao = 0. (5.7) 
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We calculate each term above separately. A direct calculation gives the first 
integral 



f „ /6mx ■ p 9m \ , 



'Si(O) 

From the area formula and Lemma |5. 11 we have 

,,a „a 



Jsiio) Jsr(p) 
h-a[ - P'')G(y) d-0 = ^^^^ + 0{R-^^n. 

JSn(v) ^ 



Moreover, the error term can be bounded by 



Since the ADM mass m > 0, we can choose p 

p{R,ij) = Ccs + e{R,ij) (5.8) 
so that the identity ()5.7p holds, where 

\e{RM < |(i + \p\) + + mm + + mntpwD^iJi). 

2. Find the solution ip by iteration. We consider the isomorphism 
L : (KerL)^ ^ Range(L) for (KerL)^ C C2'°(5i(0)) and Range(L) C 
C*''"(5i(0)). If we denote the right hand side of (|5.6p by /{pjRjip), i.e. 

f{p, R,^) = - - 1^ - R'^^G o F,^n{x) - R'+^E, o F,M^), 

we know f{p{R, ip), R, ip) € Range(L) for any R, ip with ||^||(^2,a < 1. There- 
fore, any ipQ with H^qIIc^." ^ 1; there is ipi G (KerL)-*- such that 

L^p^ = f{p{R,^Po),R,^Po). 
Moreover, we use the Schauder estimate and the fact that ipi £ (KerL)^, 
||V'l|lc2>" < c||/(p(i?,V'o),^,V'o)||co,c« 

< ^i\Ccs\ + 1) + ^^dlV'ollc^,'^ + ||^0||?,2,.) 

< ^(1^-1+1)+ f- 

For any R large enough (independent of ipo), we have || "0111(^2,0 < 1. We 
continue the iteration process and get a sequence of functions {ipk}'kLo 
isfying 

Ltpk+i = f{p{R,'ipk),R,Tpk) and ||V'fc+i||c2>- < 1- 
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The Arzela-Ascoli theorem says that there exists -000 £ C'^'^{Si{0)) such 
that a subsequence of {tpk} converges to tpoo hi C^'^ for < ^ < a. Moreover 
ipoo is a solution to (|5.6p . 

L^poo = f{p{R,Tpoo),R,1poo)- 

Let (j)oo{y) = "^00° F~^{y), then the surface Mr = {z : z = y + R~°-(l)oo'^,y G 
Snip)} has constant mean curvature equal to (2/i?) — {Am/E?). 

To complete the proof of Theorem [2l we need to compute 

lim ^ ■ (5.9) 

By the uniqueness of the constant mean curvature foliation, {Mr} are equal 
to those constructed in [HY96j , and therefore ()5.9p converges to the Huisken- 
Yau center of mass C^y. We now prove that (j5.9p also converges to the 
Corvino-Schoen center of mass C^g. Let F be the diffeomorphism defined 
by F{y) = y + R~"-(j)ooi^, then 

where JF is the Jacobian from the area formula, JF = 1 + 0(i?~^~°). Now 
we can use the fact that the area of Euclidean sphere is 0{R^) and the 
estimate for the center p in (|5.8p to conclude 

= Cg5 + e(i?,(/>oo) + 0(i?""). 

Therefore, after taking limits, the Huisken-Yau center of mass C^y is equal 
to the Corvino-Schoen center of mass and, therefore, is equal to Cf. □ 

Acknowledgments 

I would like to thank my thesis advisor Prof. Richard Schoen for suggest- 
ing this problem and for the remarkable ideas and comments he provided. 
I also would like to thank Simon Brendle and Damin Wu for very useful 
discussions and Justin Corvino for his interest in this work. 



References 



[B86] Bartnik, R., Mass of an asymptotically flat manifold, Comm. Pure and Appl. Math, 
Volume 39 (1986), pp. 661-693. 



32 



LAN-HSUAN HUANG 



[BLP03] Baskaran, D., Lau, S. R., and Pctrov, A. N., Center of mass integral in canonical 
general relativity. Ann. Physics 307 (2003), no. 1, 90-131. 

[B087] Beig, R. and O Murchadha, N., The Poincare group as the symmetry group of 
canonical general relativity. K\m. Physics 174 (1987), no. 2, 463-498. 

[CS06] Corvino, J. and Schoen, R., On the asymptotics for the vacuum Einstein constraint 
equations, J. Differential Geom. Volume 73, Number 2 (2006), pp. 185-217. 

[CW08] Corvino, J. and Wu, H., On the center of mass of isolated systems, Class. Quan- 
tum Grav. 25 (2008). 

[HY96] Huisken, G. and Yau, S.-T., Definition of center of mass for isolated physical 
systems and unique foliations by stable spheres with constant mean curvature, Invent. 
Math, Volume 124 (1996), pp. 281-331. 

[RT74] Regge, T. and Teitelboim, C, Role of surface integrals in the Hamiltonian for- 
mulation of general relativity, Ann. Phys. Volume 88 (1974 ), pp. 286-318. 

[SY79] Schoen, R. and Yau, S.T., On the proof of the positive mass conjecture in General 
Relativity, Corrun. Math. Phys. Volume 65 (1979), pp. 45-76. 

[SY81] Schoen, R, and Yau, S.T., Proof of the positive mass theorem II , Comm. Math. 
Phys. Volume 79 (1981), pp. 231-260. 

[Y96] Ye, R., Foliation by constant mean curvature spheres on asymptotically at manifolds. 
Geometric Analysis and the Calculus of Variations, Int. Press, Cambridge, MA, (1996), 
pp. 369-383. 



